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ASYMPTOTICS OF THE PARTITION FUNCTION FOR RANDOM MATRICES VIA 
RIEMANN-HILBERT TECHNIQUES, AND APPLICATIONS TO GRAPHICAL 

ENUMERATION 

N. M. ERCOLANI AND K. D. T-R MCLAUGHLIN 


Abstract. We study the partition function from random matrix theory using a well known connection to 
orthogonal polynomials, and a recently developed Riemann-Hilbert approach to the computation of detailed 
asymptotics for these orthogonal polynomials. We obtain the first proof of a complete large N expansion 
for the partition function, for a general class of probability measures on matrices, originally conjectured 
by Bessis, Itzykson, and Zuber. We prove that the coefficients in the asymptotic expansion are analytic 
functions of parameters in the original probability measure, and that they are generating functions for the 
enumeration of labelled maps according to genus and valence. Central to the analysis is a large N expansion 
for the mean density of eigenvalues, uniformly valid on the entire real axis. 
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1. Motivation and Background 

In this paper we will consider asymptotics for the following family of integrals: 


( 1 . 1 ) 


( 1 . 2 ) 


(^1 : ^2 : ■ ■ • 5 ^ 1 / ) — 




-N^ 

JR 

JR 

[ [ 


N 


i=i 


d^X, 


y(A; h,...,U) = Ft (A) = ViX) = Ja^ + ^ 4A^ 


fe=i 


where the parameters {ti ,..., t^} are assumed to be such that the integral converges. For example, one may 
suppose that u is even, and > 0. We will use the following set for allowable t = {ti,... ,t^). For any given 
T > 0 and 7 > 0, define 


T(r, 7 ) = {t€R^: ltl<T, 

j=i 


There are three main results in this paper. Theorems |1.1| , |1.2| , and |1.3| . The first of these is the following 
theorem. 


Theorem 1.1. There isT > 0 and 7 > 0 so that fort € T(T, 7 ), one has the N —> 00 asymptotic expansion 
(1-3) log = ^^eo(t) + ei(t) + ^£ 2 ( 1 ) + • • • . 

The meaning of this expansion is: if you keep terms up to order N~‘^^, the error term is bounded by , 

where the constant C is independent oft for all t G T(T, 7 ). For each j, the function 6 ^( 1 ) is an analytic 
function of the (complex) vector t, in a neighborhood ofO. Moreover, the asymptotic expansion of derivatives 
o/log(ZAr) may be calculated via term-by-term differentiation of the above series. 

Remark: In the statement of the Theorem, t S T(T, 7 ). This is not the largest domain where the 
asymptotic expansion holds true. What is really required is the existence of a path F in connecting t 
to 0 such that for all t G F, the associated equilibrium measure (see Section ||) is supported on a single 
interval, with strict variational inequality off the support, strict positivity on the interval of the support, and 
square-root vanishing at the endpoints. A global characterization of a maximal domain where the expansion 
holds true would lead us too far from the main focus of this paper, and we will not pursue this here. 

The integral Z^fti, ■ ■ ■ ,t,^) = Zi\f{t) appears in a number of areas of mathematics and mathematical 
physics. We will explain some of the motivations for studying the asymptotics of Z^lt) in the following 
subsections. In what follows it will be useful to define the ratio: 


(1.4) 


^Ar(t) = ZAr(t)/ZAr(0). 


1.1. Motivation: statistical mechanics of a log-gas. The integral 0) has the natural interpretation 
as the partition of function for a statistical mechanical system of particles on the line, with logarithmic 
interaction potential, in the presence of an external field whose potential is Vt (A). The asymptotic behavior 
for N 00 may then be interpreted as the limiting behavior of the statistical mechanical system in the low 
temperature and many particle asymptotic limit. 

With N fixed, the integral Zj^ifi,... ,t,f), and indeed all relevant statistical observables, may be expressed 
in terms of families (parameterized by {ti,..., t,ff, and N) of associated orthogonal polynomials (see Section 
1.4 below). By itself, this observation is not terribly useful. However, it turns out that these families of 


orthogonal polynomials are in a certain sense completely integrable. Indeed, all limiting asymptotic questions 
about the orthogonal polynomials (and hence concerning the statistical mechanical system) are in principle 
explicitly computable (see Section ^ below). 
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1.2. Motivation: random matrix theory. The integral ( |l.l[ ) is also of fundamental importance in the 
theory of random matrices. As is well known (see, for example, the review text [^ ), in the theory of random 
Hermitian matrices from the so-called ” Unitary ensemble”, one considers the measure on N x N Hermitian 
matrices given by 


(1.5) 


d^t = — exp{-A^ Tr [Vt{M)]}dM, 
Zn 


where dM is Lebesgue measure on the matrix entries, i.e. 


N 


dM = Y[ dMj 

j<k j=l 


33 ’ 


R T 

where denotes the real part of the matrix entry Mjfc, and denotes the imaginary component of the 

matrix entry It is a basic fact that the measure (1.5) induces a probability measure on the eigenvalues, 
with density 


( 1 . 6 ) 


Zn 


exp -N^ 


N 


^ Vt(Aj) + \Xj - Xi\ 


N m 

i=i 3¥^i 


d^X, 


where Zn = -^Ar(t) appeari ng i n (|l.6| ) is precisely the partition function defined in (IT). Thus asymptotics 
such as those contained in (y_j) yield asymptotic information concerning the statistics of the eigenvalues of 
these random matrices. For example, by differentiating log^Ar, one obtains 

(1.7) ^logZv=-A^E(TrM^), 

where E denotes the expectation with respect to the probability measure dfit- Now if one has established 
O’ then in conjunction with ( |1.7| ), one learns the following: 

(1.8) lim e(^Tt M^) =^eo{h,...,U). 

Similar calculations show that o implies 

(1.9) lim {E ( Tr (M") . Tr (M”*)) - E ( Tr M"*) . E ( Tr M")} 

Y—too 


52 


dtmdt. 


-eo{ti, ... ,tv). 


Observe that (1.8) and Q together imply that the fundamental random variables {;^Tr are 

asymptotically uncorrelated. 

1.3. Motivation: Graphical Enumeration. 

1.3.1. The Viewpoint of Gaussian Expectations. One way to try to think about the partition function Zn 
in (1.5) is as a Gaussian expectation of the interaction term: exp{—iV Tr [Vt{M) — This is the 

viewpoint taken in much of the physical literature on random matrix theory, for example . 

To better explain this we briefly review some facts about Gaussian measures: 

A measure ^ on R” is called Gaussian if its characteristic function, ^(k) := has the form 


( 1 . 10 ) 


^(k) = exp{i(p,k) - -(Qk,k)}, 


where p is the mean and Q is called the covariance of p. For mean zero and nonsingular covariance, 


( 1 . 11 ) 


dp = exp{ —-(Ax, x)}dx 


where A = Q~^. We consider dp of the form (1.11). 

The salient feature of Gaussian measures, and the one which is fundamental for applications to graphical 
enumeration, is that expectations of general polynomial functions can be reduced to quadratic expectations. 
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The explicit recipe for this, generally referred to as the Wick formula, states that for linear functions li on 

R’", 


( 1 . 12 ) 


W2...4fe) 




r2^S2. 


■ ■ ■ 


where (•) denotes expectation with respect to d/r and the sum is taken over the {2k — 1)!! (= {2k — l){2k — 3) 
... 1) Wick couplings of 1, 2,..., 2k. A Wick coupling is a partition of 1, 2,..., 2fc into couples {ri,Si) such 
that ri < r 2 ■ ■ ■ < Tk and Si > Vi. The quadratic expectations are completely determined by {xipj) = qij 
where qij is the ij*^ entry of the covariance matrix Q. Expectations of odd polynomials vanish. 

The proof of the Wick formula follows from a comparison of the Taylor coefficients of the characteristic 
function ())(k) and its logarithm. We refer the reader to |23, P. 9] for details. 

The random matrix measure dehning the Gaussian Unitary Ensemble (GUE) and given by 


(1.13) d/x = 2-^71"^ exp j-iTrM^jdM 

is manifestly a Gaussian measure. It is straightforward to work out the covariance matrix from which one 
may conclude that 

(1.14) {mijmj^) = 1, and (wj^mfe;) = 0 for {i,j)^{k,l). 

1.3.2. Enumerating Maps. There are by now a number of striking examples which illustrate the power of 
random matrix methods for calculating explicit solutions of combinatorial problems related to graphical 
enumeration. One of these is the problem of enumerating maps. 

A map is a graph which is embedded into a Riemann surface so that 

(1) the (images of the) edges do not intersect; 

(2) dissecting the surface along the edges decomposes it into a union of open cells; these cells are called 
the faces of the map. 

Given this definition, one can give a precise combinatorial description of a class of maps in terms of edge 
identifications between a collection of faces. For instance one can pose the question of how many maps can 
be constructed from a single face, having 2k edges around its boundary, by identifying its edges in pairs. 
It is straightforward to see that this number is {2k — 1)!!. A more subtle question is to ask how many of 
these maps lie on a surface of genus g. It turns out that this counting problem has a direct and natural 
reinterpretation in terms of the combinatorics of Wick couplings. Let Sg{k) denote the number of one face 
maps with 2k edges on a Riemann surface of genus g. A generating function for these numbers can be 
directly expressed in terms of a random matrix moment fl^: 


[k/2] 

(1.15) (TrM^'^) = ^ eg{k) 

9=0 

There are extensions of these kinds of calculations to more general classes of maps along with remarkable 
applications to the calculation of geometric invariants of moduli spaces of Riemann surfaces. We refer the 
reader to for a good general description of these results. 

1.3.3. Diagrammatic Expansions. We return now to considering the evaluation of the partition function for 
the deformed unitary ensemble Zjv(t). Here we will summarize some beautiful work relating the asymptotic 
expansion (0 to some problems in enumerative geometry. As mentioned above, in the authors asserted 
the existence of the expansion (|l.3| ), and presented a very elegant consistency argument, which is described 
in subsection |l.5[ The components of t are viewed as deformation parameters; when t = 0 one has the 
partition function of the original Gaussian Unitary Ensemble. 

To fix ideas, will will restrict our attention to the case considered in |3j . So in this subsection, we will set 
all parameters tj = 0, except for t = t^. (While in ||] the authors considered more general deformations, the 
authors gave a more detailed discussion for this case.) 
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Consider the Taylor series expansion corresponding to all parameters tj = 0 except for t = 


(1.16) 


exp 


N 


= yl(zl 

n ) \ N 


n>0 


{TrUy 


This is clearl y globally convergent for all t. The difficulty arises when one considers the Gaussian expectation 
(using (1.13)) of this exponential function. Commuting the expectation integrals with the sum on the right 
hand side, |^] produced the formal Taylor series expansion of the partition function around t = 0: 


(1.17) 


ZN(t) = (exp<{-—TrM 


Zn{0) 


Zjv(0,0,0,t4 = C0,...,0) 


n>0 ■ ^ 


The quotation marks indicate where one has proceeded formally. Indeed it is manifest that this is not 
a convergent series expansion since the integral corresponding to (exp{ —(t/iV)TrM‘^}) converges only for t 
with positive real part. Properly speaking, this is really shorthand for a series of identities relating derivatives 
oi Zm evaluated at t = 0 to expectations of powers of Tr M'^. 

Again, use of the Wick calculus allows one to replace the matrix moments appearing in the above formal 
series by a generating function for enumerating a certain class of labelled 4-valent maps (or disjoint unions 
of maps) which we will refer to as diagrams. A four valent diagram consists of 

(1) n 4-valent vertices; 

(2) a labelling of the vertices by the numbers l,2,...,n; 

(3) a labelling of the edges incident to vertex a (for a = 1,..., n) by letters v, jtr, ka, ia (This alphabetic 
order corresponds to the cyclic order of the edges around the vertex); 

(4) a partitioning of the labelled edges into pairs. 

(Once the number n is specified, the set of all n-vertex, four valent diagrams is in one-to-one correspondence 
with the set of all pairings.) Connecting the edges according to the pairing yields a graph, together with a 
cyclic ordering of the edges around each vertex. For definiteness, we will adhere to the rule that this cyclic 
ordering corresponds to a clockwise orientation of the edges around a vertex. 

If the underlying graph is connected, the cyclic ordering determines a map associated to the diagram. In 
that case we will call the diagram a g-map if the map is on a Riemann surface of genus g. (Note that a 
g-map as defined here has more structure than a map as defined in the previous subsection since a g-map 
carries a labelling of its vertices and edges.) If the graph is not connected, then the edge labelling associates 
a map to each connected component, and the conglomerate of maps will be called a g-diagram. 

The use of the Wick calculus yields the following different representation of (1.17), which clearly demon¬ 
strates that the partition function is a generating function for the enumeration of diagrams: 

1 /-t" 


(1.18) 


^iv(t)“ = ” 


V- 

^ n\ 

n>0 


N 


ff{A-vdlent, n-vertex, g-diagrams}A^^ 2 g+n^ 


Exploiting the relationship between the terms of a Taylor series and those of its logarithm one can write 
down an equally formal representation for the the logarithm of the partition function which can be regarded 
as a generating function for connected n-vertex diagrams: 


(1.19) 


logiAr(t)“ 


1 / _ t\^ _ 

” ~! ( ^ #{4-valent, n-vertex g-maps}iV^“^®^". 

n>o ^ g>0 


Finally one can make another leap beyond rigor and resum the terms of the previous formal series, as was 
done in 1^, to formally order the series by the genus of the surface to which the diagram maps. This yields 
what is referred to in the physics literature as the genus expansion, 


logZjv(t)“ = ” yEg{f)N^-^\ 

9 


( 1 . 20 ) 
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where Eg{t) = Un>i is a formal series (possibly convergent) in which each of the coefficients 

Kg{n) is the number of connected maps of genus g with n vertices (all 4-valent). 

The authors of fully appreciated that this intuitive “deduction” was completely formal and in fact 
suggested an approach for providing a rigorous derivation. This approach was never fully pursued (we 
briefly review its history in Section 1.5). 

We have taken a differe nt ap proach and our proof of Theorem 1.1 provides a rigorous validation of the 
form of the representation ( 1.20 ). When restricted to the four-valent case, our result may be summarized as 
follows. 

A. There is a positive number T 4 such that for 0 < t < T 4 , the partition function ZMit) possesses the 
following asymptotic expansion. 

(1.21) \ogZM{t) = eo(t) -I- N~'^ei{t) + N~‘^e2{t) H-. 

The meaning of this expansion is: if you keep terms up to order N~'^^, the error term is bounded by 
The coefficients eg{t) possess analytic continuations to a neighborhood of t = 0. 

B. The coefficients eg(t) are related to the counting of 4-valent g-maps via the following formula: 

9" 

(-l)”^eg(0) = #{4-valent, n-vertex, g-maps}. 

For the case of more general deformations, and the full asymptotic expansion (Ol), we have the following 
theorem concerning the coefficients eg{ti .. .tn). 


Theorem 1.2. The coejficients in the asymptotic expansion (l.c) satisfy the following relations. Let g be a 
nonnegative integer. Then 

1 


( 1 . 22 ) 


i{ti .. .t^) — 


n,- > 1 


m!... nj. 


(-ti)"i ... (-ti.)”‘'«:g(ni, ...,n^) 


in which each of the coejficients Kg{ni,... ,n,^) is the number of g-maps with nj j-valent vertices for j = 

1.4. Motivation: The Connection to Orthogonal Polynomials. As mentioned above, the integral 
(0 is intimately connected to the theory of orthogonal polynomials. Consider the measure 

(1.23) W]s[{x)dx := exp [—NVt{x)]dx. 


Let us define {pj{x; N, t)}^g to be the sequence of polynomials orthogonal with respect to the measure 
WN{x)dx. That is, {pj{x-, N, t)}^Q satisfies 


(1.24) 


PjPkWNdx = 


0 j 

1 j = k ’ 


and Pj{x; N, t) = + •••, > 0. (The leading coefficient 7 ^"^ is of course dependent on the 

parameters ti,...,tv\ however, we suppress this dependence for notational convenience.) The fact of the 
matter is that ^Ar(t) may also be defined via 


(1.25) 


N-l 

Z^(t) =m Y[ (yf 

e=o 


Zn is also defined via 


Co 

Cl • • 

CN-1 

Cl 

C2 • • 

Cat 


(1.26) 


ZN{tl, . . . , tjy) = A^! 


cat-i Cat 


C 2 N -2 
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where Cj = J^x^WNix)dx are the moments of the measure WN{x)dx, and the determinant ab ove is called a 
Hankel determinant (see, for example, Szego’s classic text [|^). The asymptotic expansion (0) constitutes a 
version of the strong Szego limit theorem for Hankel determinants. The strong Szego limit theorem concerns 
the asymptotic behavior of Toeplitz determinants associated to a given measure on the interval (0, 27r). We 
refer the interested reader to 1^ for more information. 


It is well known in both the approximation theory literature and the random matrix theory literature 
(see, for example, and the references therein, or and references) that ^^^(t) satisfies the following 

leading order asymptotic behavior: 

J™ . .t„) 

N^oo ^ 


sup|- [ Vt{X)dfj.{X) + [ [ log\X - r]\dfi{X)dfj.{r]) 
fj.GA L J J J 


where A is the set of all positive Borel measures on the real axis, with unit mass. The variational problem 
posed above is referred to in the approximation theory literature as the problem of determining the equilib¬ 
rium measure for logarithmic potentials in the presence of an external field (see ||^). It is well known that 
the supremum is achieved at a unique measure fi* (see, for example, p^), and that for real analytic external 
fields, the equilibrium measure ji* is supported on finitely many disjoint intervals, and on the interior of each 
interval, the equilibrium measure has analytic density |^. 

1.5. Brief discussion of history. In |^, Bessis, Itzykson and Zuber considered asymptotics for the integral 
(0), but for the special case in which tj = 0 for j ^ 4. That is, they considered the integral 


(1.27) 


zT\t) = 


exp -N^ 




N E 

j=i ^ ^ 


d^X. 


They computed asymptotics for this integral by using (1.25), together with some reasonable assumptions on 


the asymptotics of ratios of the leading coefficients yj” . The basic idea is to re-write ( 1.251 ) in yet another 
form, namely 

(1.28) 


^^ogZr\t) = 


^ log (iV!) - ^ log + E ^^^2 log (^i-O ’ 

1=1 


where denotes the recurrence coefficient associated to the orthogonal polynomial sequence, xpj{x) = 
bjPj^i{x) + ajPj{x) + bj-ipj-i{x). It easy easy to see that the recurrence coefficients are related directly to 
the leading coefficients y]^\ via bj = yj^V7j+i- 

The authors then made the following asymptotic assumptions concerning the recurrence coefficients : 


(1.29) 


(feW) = ro ( ) -h —r2 ( 


iV’ 


1 

N2’ 


l_ 

TV’ 


7V2fe"’2fc 1 +••■ 


where the functions r 2 i{x,t) are assumed to be infinitely differentiable functions of the variable x. With 
these assumptions in hand, Bessis, Itzykson, and Zuber then substitute ( |1.2g| ) into ( |1.28D , and apply the 
Euler-Maclaurin summation formula to replace each summation by a series of integrals. 

Following that, there is an argument that if an asymptotic description such as (1.2!\) holds, then one may 


deduce the form of the functions r 2 i. This is achieved by studying a nonlinear recursion relation satisfied 
by the recurrence coefficients b^f^\ This in turn leads to a prediction for the asymptotics of The 

authors compute explicitly the first three terms, and present a conjecture for all terms. Their conjecture is 
that for j > 2, e^^^'^\t) is a rational function of t, with one pole in the finite plane, of order 5(j — 1), located 
at t = —1/48. 
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Quite separately, a rigorous analysis of the recurrence coefficients for polynomials orthogonal with respect 
to weights of the form dx [ p0[ was carried out by analyzing the very same nonlinear recursion relation. 

This yielded a proof of the leading order information contained in ( 1.29| ), but valid only if j ^ oo. 

On the other hand, in |^, the authors use (1.29) as a uniform expansion, valid even for j = 1, as fV ^ oo. 
It would seem plausible that such an expansion would hold for j oo, N oo, such that j/N ^ a: > 0. 


However, it is entirely possible that for j such that j/N approaches 0, there is some correction to (1.29) 
which may affect the asymptotic computation of at some finite algebraic order in N. 

For this reason, it is of great interest to find some independent way to compute the asymptotics for 
log(Z.(t„...,t.)). ^ 

In a different direction, in Q, techniques from the theory of isomonodromic deformations, as well as tech¬ 
niques for the asymptotic analysis of Riemann-Hilbert problems were applied to the problem of determining 
asymptotics for the polynomials Pniz] N, t) orthogonal with respect to the 1-parameter family of measures 
exp [—N{x'^/2 + tx'^)]. The results were used to prove that the associated random matrix model obeys the 
universality conjecture. That is, it was proved that the local statistics of the eigenvalues of random matrices 
with the probability measure exp [—iV Tr {M'^/2 + tM'^)]dM converge as ^ oo to a universal random 
point process, whose correlation functions are given in terms of the so-called sine kernel. 

In addition, in recent techniques for the asymptotic analysis of Riemann-Hilbert problems were 

applied to the problem of determining asymptotics for a wide class of orthogonal polynomials, which contain 
the orthogonal polynomials Pn{z', N, t) as a subset. The results were used to establish this universality 
conjecture of random matrix theory, for a family of real analytic probability measures. As a by-product of 
the analysis in asymptotics for the recurrence coefficients were obtained. 

The upshot of this is that there is emerging a new way to compute asymptotics of log (Zjv(t)). In both Q 
and [^-[11, an important asymptotic technique is the Deift-Zhou steepest descent / stationary phase method 
for the asymptotic analysis of Riemann-Hilbert problems. This method was introduced in |n], and further 
developed in and |p^ . 

One way to proceed is to compute the asymptotics for the orthogonal polynomials, and then deduce the 
asymptotics for the recurrence coefficients However, this still requires an analysis of the recurrence 

coefficients which is uniform in j, i.e. for j = 0 through j = N — 1. This appears to be an onerous task, and 
so we shall not adopt this approach here. Instead, we will proceed, starting in the next subsection, by the 
direct analysis of yet another representation of log (.^jv) (formula O). 


1.6. A statistical mechanical formula, and orthogonal polynomials. In this paper we will analyze 
the large N behavior of Z]\{{t) by direct analysis of the formula (1.33) for the logarithmic derivative of 
The derivation of this formula proceeds as follows. One begins by computing the logarithmic derivative of 
Zn- 


(1.30) 





N 




1=1 


X exp < — 


N 


d^X. 


^ ^ Vt{Xj) - ^ ^ log |A, - A^l 
i=i 

Of course, this has the immediate interpretation as an expectation value, with respect to the probability 
measure (0: 

8 ( ^ 

(1.31) — log(Zjv) =E^ 

Now it is well known (see, for example, |2^) that expectation values of functions of the XjS can be expressed 
in terms of the associated orthogonal polynomials. In particular, the mean density of the Aj’s is given by 
the so-called ” one-point function”, 

1 


Pa (A) = — exp[-At4(A)] ^ Pi(A; N,tf. 


(1.32) 
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Using this connection, one has the following remarkable formula for the partial derivatives of log (Z]\f): 

(1.33) ^ log (Zn) = (^1 Tr = -N^ J X^p^^\x)dX 

in which the first equality follows from (O)- The fundamental theorem of calculus readily implies 

(1.34) Ziv(t) = Zjv(0)exp|-Af2 / [ p[^\x)^tVdX ■ m] 


Since Z^it = 0) is explicitly known {Zpf(0) = •••), the formula (1.34) shows that if we have global 
asymptotics for p\^\ we can deduce an asymptotic expansion for Zpf{t). We will analyze the one-point 
function, globally on the real axis, and determine a uniform asymptotic representation for it. Then we will 
evaluate the integral appearing in the right hand side of ( [l.33| ), and from this obtain a uniform asymptotic 
expansion for (1.34). This is the procedure by which we shall compute the asymptotic expansion of logZf^. 
The following theorem is a fundamental result which, following the above prescription, yields Theorem 

uTI. 


Theorem 1.3. There is T > 0 and 7 > 0 so that for all t G T(T', 7 ), the following expansion holds true: 

/ OO 

/(A)pW(A)dA = /o + iV-Vi + N-^f2 + ■■■ , 

-00 

provided the function /(A) is C°° smooth, and grows no faster than a polynomial for X —> 00 . The coefficients 
fj depend analytically on t fort G T(T, 7 ), and the asymptotic expansion may be differentiated term by term. 

In the statement of the above Theorem, the meaning of the asymptotic expansion is this: if you stop after 
the term N~‘^^fj, the error term is 

Remark 1.4. We note that in some very special instances results of the above form may be deduced from 
known results in the literature. For instance, taking the Gaussian limit of (1.33) one has 

r) 1 

■log(Z^)|t=o = -iV2 / A2V^^(A,t = 0)dA = -iV2(- 

J —00 * 


dt2k 


AT-fc Tr 


The last term, —Tr M^^), which by formula ( 1.15 ) is equal to a finite series, —Sgio ^gi^) {^) > 
is a genera ting function for counting one-face maps with k edges. This series is manifestly even in N. In the 
s ettin g of (1.33| ) this series is actually a generating function for counting 2fc-valent one-vertex maps (see also 
(^- (^. 261 ) ). There is no inconsistency in this since by duality these two countings are equivalent. 

Rema rk 1.5. Additionally, in Q, Albeverio et al have recently considered asymptotic expansions of the form 

( 1.35 ), for the special function / 2 (A) = (z — X)~^, under the assumption that the external field is even, but 
they allow for the equilibrium measure to be supported on one or two intervals. They assert an asymptotic 
expansion for such integrals which is in inverse powers of IV; our analysis shows that there is a complete 
expansion, in inverse powers of JV^. 

The paper is organized as follows. In Section ^ we present a self-contained description of the connection 
between random matrix integrals and the enumeration of maps on Riemann surfaces, according to genus and 
valence; we also prove Theorem 1.2, using Theorem 0 In Section ^ we summarize the procedure developed 
in |^-|^ to establish rigorous uniform global asymptotics for the associated orthogonal polynomials. In 
Section 1^ we use the results of Section || to derive exact formulae for pj^^(A), the mean density of eigenvalues 
at finite JV, in terms of the Riemann-Hilbert procedure. Explicit in the formulae derived is a complete 
asymptotic expansion for . In Section we use the formulae derived in Section ^ to establish Theorem 

01. 


2. Enumeration of g-maps 

In this section we will explain the connection between the evaluation of the Gaussian moments ( 0^=1 
(TrMl)'^^), and g-diagrams and g-maps (as defined in subsection 1.3.3| ). Following that we will prove The¬ 
orem 1.2 which states that the asymptotic expansion ( |I.3| ) is a generating function for g-maps, enumerating 
them according to the number of vertices with different valences and the genus. 
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2.1. Gaussian matrix integrals and enumeration of diagrams. For simplicity we first discuss in detail 
the case of the pure moment (TriW^)”; afterward we will indicate how this extends to general mixed moments. 
The expression (TrM"*)" has the form 


( 2 . 1 ) 






where the sum is taken over all configurations, where by a configuration one means an assignment to each 
of the 4n labels appearing in the list {v, ju,kv, a value from {1,..., N}. There are thus config¬ 

urations. Each configuration corresponds to a selection of 4n matrix entries. 


We compute the expectation of this quantity with respect to the original probability measure (1.13), 
interchanging integration with summation: 

( 2 . 2 ) 

— ^ ^ ^ki ,j2 ^j2,k2 ■ ■ ■ ,jn ^jn,kn ,^n ) 

configs 

By the Wick calculus, the expectation of each te rm in the above sum can itself be evaluated as a sum over 
Wick couplings as they were described in Section 1.3.1 . To appreciate this one should first observe that each 
of the matrix entries appearing in the expectation on the right hand side of (^^) is a linear function. Next, 
we introduce the following notation: 


(2.3) 

where 

(2.4) 

and for u = 2, 

(2.5) 

( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 




1 .^1 ^^1 j'^l .U ^*2 02 ,^2 ■ ■ ■ On ) 


,n. 


= (/ 1/2 ■ • • fin) , 

fl = /2 = mj-^^kn fs = fi = , 


f4(v — l) + l = ''TT-inOn) 
/4(l^-l)-|-2 = 
/4(y-l)-|-3 = 
f4(v — l)+4 ~ 


The definitions ( 2.4)-( p.^ ) implicitly define an invertible mapping (r, s)(-) = (?■(•), s(-)) between the inte¬ 
gers {1,... 4n} and the matrix indices {(v, jf/), [jv, ku), (k^, (f^, i,^)} which we will describe as follows: 

(2.9) (r,s)(-) = (r(-),s(-)) : {1, ■ • ■, 4n} ^ {(v, >), (>,/c^,), (A;,,, C), (4, v)}"^i, 

k ^ (r(fc),s(fc)), so that fk = mr(k),s{k) (using (^-(^). 

We now apply the Wick calculus: 

(2.10) ^ (/l---/4n) 

configs 

= E E /a;i,2 ) {/^2, 1 /^ 2 , 2 ) * ’ * {/t.^2n,l/^2n,2 ) 

configs “£W4n 

In this setting a Wick coupling, lo = ujij, i = 1,..., 2n, j = 1, 2, is a partition of 1, 2,..., 4n into 2n couples 
(wip, Wi_ 2 ) so that uji^i < W 2 ,i < • • • < W 2 n,i, and u}i ^2 > for i = 1,..., 2n. We denote the set of all such 

Wick couplings of the integers 1,..., 4n by IT 4 „. _ 

Now the outer summation over configurations appearing on the right hand side of ( ^.10 ) is over the 
different realizations of the variables /i,...,/ 4 n determined by the configuration. We may interchange 
orders of summation: 
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( 2 . 11 ) 


((TrM4)")= ^ (/i---/4n) 

configs 

= E E /i2J2,2 } ’ ’ ’ {/)tJ2n,l/i2J2Ti,2 } ' 

i^eW 4 ^ configs 


^1,1 f^1,2'} ^2,if<^2,2') ' ‘ ‘ {/^ 2 ti , 1 /)^ 27 i ,2 } 


The individual terms 

( 2 . 12 ) 

are either 0 or 1. For each Wick coupling oj, certain configurations contribute a 1 to the inner summation, 
and other configurations do not contribute. The question of evaluating the original integral ((Tr M^) ) has 
been reduced to computing, for each Wick order ing uj € W 4 n, the number of configurations providing unit 
contribution to the inner sum appearing in ( 2.11 ). For each different Wick coupling oj, it turns out that this 
number is of the form and accepting this for the moment, we have the formula 

(2.13) ^ (fuiijuii,,) 
configs 

To determine the contribution for each Wick coupling, and to explain the quantity F{uj), let 

us first recall the direct relationship (2.4)-(2.8) between {fj} and the original variables of integration 

; and the identities ( 1.14 ). Using the mapping (r, s)(.) defined in (2.9) 
, the inner sum in ( | 2 . 1 l| ), over all configurations, is of the form 

(2.14) ^ ^ (™r((.iJl,l),s(wi,l)™r(wi, 2 ),s(l^l, 2 )) (™r'(w 2 ,l),s(l^ 2 ,l)^l'(i.J 2 , 2 ),s(l.J 2 , 2 )) 
configs 

{l^r(0J2-n,l),s{iJ2n,l)'^r{iJ2n,2),s{uJ2rv,2)') ' 


Each different configuration appearing in (2.14) is a different choice of the matrix indices {(iv,ju), 

(t'l,,^)}, and for each one we then ask if the contribution is 0 or 1. Using ( 1.14| ), the contribution 
will be 1 if and only if for each ^ = 1 ,..., 2 n, 


(2.15) 






^) = 1 > 


and this is true if and only if 

(2.16) ?'(w/ 2 ,i) = s(a;^, 2 ) and s(w^,i) = r(a;^, 2 ) for each /r = 1 ,..., 2 n. 

The sum (2.14) is precisely the number of configurations for which the 4n equalities (H) hold true, 
and since r(a;^p) is, for each n, one of {ii/, (and similarly for s{u!fj,^ 2 ), s{uj^^i), and )> it 

follows that the sum over configurations, is now over 

jy-i — 1, . . . At, 


so that ( 2.16| ) holds. 


The calculation of (2.14) now proceeds as follows. For each Wi ck co upling uj we write out equalities 
amongst the 4n indices of summation enforced by ( 2.16| ) and group them in “chains” (for 

example, one might have 13 = js = £7 • • •). Each index of summation appears in two equalities (because 
we are computing ((Tr see (^)) and so each distinct chain is, fact, a closed cycle. For each closed 

cycle, the original summation is reduced to a summation over a single free parameter, ranging from 1 to N. 
Finally then, the summation ( |2.14 ) is where F{uj) is precisely the number of closed cycles determined 

by the Wick coupling w, and we have explained ( ^.13 ). 

Summarizing, we have shown that 


(2.17) 


<^(TrM4)"^= 


W4ti 


where F{(jj) i s the number of closed cycles determined by u via the equalities (|2.16D . Note that the right 
hand side of (2.17) is a finite sum, as W 4 „ is a finite set. 
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A key insight of the classic paper of g was to realiz e tha t cou nting the Wick couplings could be reduced 
to listing all possible diagrams (as defined in between (1.17) and (I.IS) above). To see this one begins with 
a collection of n 4-valent vertex configurations labelled as in Figure 2.18. A given coupling determines a 


unique way of connecting the outgoing and incoming rays of each of the n 4-valent vertices, to form “roads”. 
The result of such a glueing will be a diagram. Indeed, if we interpret each road in Figure 2.18| as an edge, 
and each edge is labelled by the “outgoing” side of the road, we clearly have a diagram. (The “faces” of the 
diagram correspond 1-1 to the closed index cycles of that coupling.) Conversely, given a diagram as defined, 
which includes the labellings, we can read off directly, from the way the incident edges are paired, what the 
couplings must be. This gives a bijection between diagrams and Wick couplings. 



Figure 2.18. 

The full count of non-vanishing contributions is then given by summing over diagrams with a weighting, 
^F(w), counts the number of configurations that give rise to a nonvanishing contribution in the coupling, 
Lu, associated to this diagram. Also the exponent of the weighting, F{uj), which was defined to be the number 
of closed index cycles in the coupling is now also seen to be the number of faces, F, in the associated diagram. 
The weighting factor is . If the diagram is connected (i.e., if it is a map), this can be expressed in terms 
of the genus of the Riemann surface that the diagram maps to by using Euler’s formula. Since the number 
of vertices of a diagram is n and the number of edges is 2n (= the number of couples in a Wick coupling), 
we have 2 — 2g = n — 2n + F from which we deduce that the weighting factor for a 4-valent, n-vertex, g-map 
is N'^~‘^ 3 +n_ yyg directly extend this to the case of disconnected diagrams since the Euler characteristic 
is additive with respect to disjoint unions. This generalization results in the possibility that g can become 
negative. Thus, to each diagram we can associate a unique integer g and then we will refer to that diagram 
as a ^-diagram. 

The beautiful connection between the quantity ((Tr M"^)") and the combinatorics of g-diagrams described 
in (and discussed above) is summarized with the following formula: 

(2.19) ^{4-valent, n-vertex, g-diagrams} 

g 


We remark that this is a finite sum, and the reader may easily verify that the nonzero contributions to this 
sum come from 1 — n < g < [(n -f l)/2], where \i] is the closest integer to i less than or equal to i. A 
straightforward extension of the above analysis to the case of mixed moments shows that 


(2.20) / (Tr Afl)"’' \ = ^{g-diagrams with Uj j-valent vertices; } = 1,..., v}N^ ‘^g+T, 

\i=i / 9 

(Again we observe that this is a finite sum.) 


2.2. Asymptotics of the partition function and enumeration of g-maps. We turn now to the proof 
of Theorem 1.2. The starting point is the asymptotic expansion in Theorem 1.1, rewritten in the form 
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This is, as usual, shorthand: the right hand side is an asymptotic expansion in inverse powers of N'^. 
Exponentiating both sides of this equation we deduce that Z^eT^ eo-ei asymptotic expansion in 

inverse powers of N'^. Some algebraic manipulation shows that 


( 2 . 21 ) 


where 


Zn = exp (TV^eo + ei) | 1 + X! f 


^=1 


h 


Ph= 7 .-n 


^ k\ 

/c=l £i + . 


n 


eii+i- 


.+^k=h,lj>0 j=l 


Formula (2.21) may appear to be nothing more than a simple calculational re-representation of Zj^. However 
( 2 . 21 ) explains concretely i n w hat sense Zn itself possesses an asymptotic repres entat ion. 

It follows from Theorem O that the asymptotic series on the right hand side of ( 2.21 ) can be differentiated 
term-by-term for t g T(T , 7 ), and the ensuing derivatives may be evaluated at t = 0. Having established 
this we now rewrite (^. 2 l|) as 


( 2 . 22 ) 


1 


1 


ZNit) = exp eo(t) -b ei(t) -b 1 ^ 62 ( 1 ) -b • • • -b ., 0 . 2 6 / 1 ( 1 ) -b ... 


' ' ]y2h-2' 

In light of the above discussion, it is clear what is meant by this asymptotic expansion, which is valid for t 


satisfying the regularity conditions of Theorem 1.1. 

Differentiating ( 2.22| ) nj times in each tj and taking the limit t —> 0, we may now conclude, by comparison 
with ( ^.2C| ), that 

(2-23) 




t=o 


= ( n (Tr 


\i=i 


= ^ #{g-diagrams with Uj j-valent vertices; j = 1 ,..., 

9 

dn 


= (- 1 )^ 


... dC'^ 


exp ^®eg(t) 


^3=0 


where n = ni -b • • • -b 

The string of equalities in (2.23) states that the finite sum over g in the second line possesses an asymptotic 
expansion in inverse powers of . We learn that all of the coefficients in this final asymptotic expansion 
must vanish from some point on. Much more information can be obtained in this manner from (2.23). For 
example, by examining ( 2.23| ) for the case of a single derivative, one may verify that for i even, 

d 

(2.24) eo(0) = —#{ genus 0 diagrams with 1 £-valent vertex} 


(2.25) 

(2.26) 


die 

A 

dU 

_d 

die 


ei( 0 ) = — #{ genus 1 diagrams with 1 £-valent vertex} 


—eg(O) = 0 forg>- 


£/2 if £/2 is even 

i/2 — 1 if£/ 2 is odd 


We observe that in (2.24) and ( ^.25 ), the diagrams have a single vertex, and so we may say that (2.24) and 
(2.25) enumerate 1-vertex g-maps. 

Proceeding further, one may examine ( ^.23 ) for the case of higher order derivatives. Tedious calculations 
may be carried out which show that the Cg ’s are generating functions for g-maps (they enumerate those dia¬ 
grams that are connected). However, this approach is quite complicated, and we adopt a different approach, 
which is to use ( ^.23 ) to discuss the relation between the logarithm of the partition function and connected 
diagrams, or g-maps, that was alluded to in (1.19|). For the sake of clarity we will present this just in the 
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case of t = ( 0 , 0 , 0 , ^4 = 0 ,..., 0 ); it will be evident from the discussion that the extension to more general 

valences is straightforward. 

We observe that, from our earlier discussion of the relation between diagrams and Wick orderings, that: 


(2.27) 


g 


# {4-valent, n-vertex, g-diagrams}A^^ ^ ~ Q^[N). 


Qn{N) is a Laurent polynomial which we will refer to as the diagram polynomial for 4-valent, n-vertex 
diagrams. We introduce an exponential generating function ||l^ for diagram polynomials: 


(2.28) 




k=l 


kl 


This is a formal generating function in the same sense as was discussed for (1.18). It is by no means 
a convergent series, but rather is interpreted in the sense that one is allowed to evaluate finitely many 
derivatives of both sides at t = 0. (Below we will indicate how to remove the formality of these arguments by 
truncating the series.) We also introduce an exponential generating function for connected 4-valent, diagrams 
(i.e. maps): 


Cit) := 


fc=i 


where 

(2.29) 


Pn{N) = Y. 


NX{^ 




and Wl®"” is defined to be the subset of W 4 „ consisting of Wick couplings that correspond to connected 
diagrams. We observe that 


J.ib 

cit)cit) = E E 


n—2 k-\-^—n 
00 n-l 

= EuS 

n=2 


Pk{N)P^_k{N). 


This series is a generating function for diagram polynomials of ordered pairs (lui,uj 2 ') of connected Wick 


couplings. There is a geometrical explanation for the appearance of the combinatorial coefficient 


this expression. It is a consequence of the fact that the labels 1,.. .n get distributed over uji U UJ 2 and if 


Iwil = k there are 


different ways of doing this. 


Next, one may prove that C{t)C(t)/2 is the exponential generating function for 4-valent diagrams which 
consist of exactly two connected components. In other words, one may define an exponential generating 
function for 4-valent diagrams consisting of exactly two connected components, and prove that this generating 
function coincides, term-by-term, with C(t)C(t)/2. Subsequently, and in the same way, one may prove that 
C‘^{t)/d\ is the exponential generating function for 4-valent, n-vertex diagrams with exactly d connected 
components. One may finally conclude that, in the sense of formal generating functions, 


(2.30) 


= ElTiC^^W =exp(C(t))-l. 


d=l 


dl 


To get around the formality of these generating functions, we repeat most of the argument from ( p.2^ ) to 
(2.30), truncating the infinite series involved. So we define 

m k ^ j.k 


k^l 




(2.31) 
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with Qk and Pk defined in ( ^.27 ) and ( ^.29| ), respectively. Next we have that 
(2.32) (^) + . 


fc=2 


where the sum on the right hand side is the truncated exponential generating function for 4-valent diagrams 

(2) 

consisting of exactly two connected components (i.e. ' (N) is the diagram polynomial for 4-valent, k-vertex 

diagrams with two connected components). Similarly, the reader may verify that 


(2.33) 


1 ™ JLfc 

CUtf = E ^ > 


d\ 


k—d 


where on the right hand side, is the diagram polynomial for 4-valent, k-vertex diagrams with d 

connected components. Putting this all together, we have shown that 

(2.34) Gm{t) = -h O . 

And since m was arbitrary, this is the proper sense in which one may interpret ( 2.301 ). 

We have established the following set of relations: 

(2.35) l + Gm{-t) = ZN{t) + 0{r+^), 

(2.36) expC™(-t) = 1 + Gmi-t) + O . 

We therefore have 

(2.37) Gm{-t)= log {ZN + 0{t^+^)). 

We may now differentiate this n times (for n < m) to obtain 


(2.38) 


dn 


■C™(-<) 


t=o 


= Pn{N) = ^logZNit) 


t=0 


So we have proven that log Zn is the generating function for connected diagrams. Now using our asymptotic 
expansion for logZjv, we have 


& 


(2.39) (-!)"_ 


9=0 


t=0 


= PuiN) 

= E] = E]-^^~^®#{4-valent, n-vertex, g-maps}. 


We have now shown that eg{t) = Eg{t) in the sense that 

Wt 


(“l)"A:7r®9(0) = #{4-valent, n-vertex, g-maps}. 


By similar reasoning one can fully establish Theorem 1.2, and we will not present these details here. 

We refer the reader to for a fuller discussion of the general cumulant relations as they pertain to 
graphical enumeration. 

3. Asymptotics for orthogonal polynomials via Riemann-Hilbert methods 

The amazing fact, due to Gaudin and Mehta [Q, is that the mean density can be expressed in closed 
form in terms of the orthogonal polynomials {pj(-; N^t )}. 


(3.1) 

N-l 

pW(A) = -e-^^WEP'=(^) 

(3.2) 

= p'nWpn- 


(N) 

Tn-i 

(N) • 

7iv 

The second equality above follows from the first because of the Christoffel-Darboux formula. (See, for 
example, formula (3.2.3) in 0.) 
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We obtain global asymptotics for by carrying out to higher order the results of @-|§. Deift, 
Kriecherbauer, McLaughlin, Venakides, and Zhou established very detailed Plancherel-Rotach type asymp¬ 
totics for polynomials orthogonal with respect to under some very general assumptions on the 

function V, including the functions V considered here. In Q, the authors demonstrated that their method 
can be used to obtain a complete asymptotic expansion for the polynomials. We have carried this procedure 
out to higher order. 


Using the global asymptotic expansion of p^^\ together with (1.34), we deduce an asymptotic expansion 
for logZjv(t). 

In this section we review the Riemann-Hilbert approach for computing and establishing a global, uniform 
asymptotic expansion for the orthogonal polynomials. In |^-|^, the authors began with a matrix valued 
function Y{z), defined by 


Y = 


:^PN 

’N 






id-/ 

(N) r r 




s — z 

^-NV(s) 


, 2^ € 


The matrix Y solves the following Riemann-Hilbert problem: 

• Y analytic in C\R, 

• Y = {I + Oi^)) 

• Y has Holder continuous boundary values Y± for z e R, 

• l+ = y_(i'=-f''),forzGR. 

The authors then make two explicit transformations: Y ^ M ^ Mi. Here we will explain these two 
transformations. 


3.1. Equilibrium measure and first transformation Y—>M. . A fundamental object in the theory of 
random matrices, as well as approximation theory, is the equilibrium measure, defined as follows: 


sup 

Borel measures /^,/.t>0,/ d^—1 


-1 V{X)dpiX) + 


log|A - p\ dp{X)dp{r]) 


There is a vast literature on this variational problem, originating in the work of Gauss (see ||^ ). Under 
general assumptions on V, the supremum is achieved at a unique measure /iy, called the equilibrium measure. 
For the external field V, considered here, it is a well-known fact (see Q, or [^) that the equilibrium measure 
is supported on finitely many intervals, with density that is analytic on the interior of each interval, behaving 
at worst like a square root at each endpoint. 

We remind the reader of the definition of the set T(T, 7 ): 


(3.3) 


T(r, 7 ) = {tG 


Itl < T, 


U > 7X! 1 ^/ 1 } 


The following theorem is perhaps well understood. A proof can be constructed by combining the results 
of Q and im . 


Theorem 3.1. There is Tg > 0 and 70 > 0 such that for all 0 < T < Tq and 7 > 70 , the following holds 
true: If t G T(T, 7 ), then 


dpv = dX, 

V’(A) = ^X(a./ 3 )(A) V(A - a){P - A) h{X), 

where h[X) is a polynomial of degree v — 2, which is strictly positive on the interval [a,/3] (recall that the 
external field V is a polynomials of degree v). The polynomial h is defined by 

Stti / y/(s - a)^J{s - (3) S- z 

where the integral is taken on a circle containing (a, (3) and z in the interior, oriented counter-clockwise. 













ASYMPTOTICS OF RMT PARTITION FUNCTION 


17 


The endpoints a and fd are determined by the equations 



V'js) 

V^(s - a)(/3 - s) 


ds = 0 


sV'js) 


ds = 2 tt. 


The endpoints a(t ) and (3{t ) are analytic functions of t which possess smooth extensions to the closure of 
T(T, 7 ), and satisfy —a(0) = /3(0) = 2. In addition, the coefficients of the polynomial h{X) are also analytic 
functions of t with smooth extensions to T(T, 7 ), with 

h{X,t = 0) = 1. 

In this paper we will assume that t G T(T, 7 ), for T small enough, and 7 large enough, so that Theorem 

7 > max{ 7 „, 7 / 3 } then t G T(T, 7 ) implies more: the neighborhoods of z = /3 and z = a defined in this 
section can be taken so large that their intersection contains a subinterval of the real axis. We will not 
present a proof of Theorem |3.1| here, but the basic idea is as follows: for t = 0, the equilibrium measure is 
explicitly known, = X(- 2 , 2 )(^)'s /4 — A^/(27r). Taking T small implies that the endpoints of the support 
cannot stray too far from ± 2 ; taking t^ larger than 7 X]j=i \^j \ implies that h should be strictly positive on 
R. 


3.1 holds true. In Section W we define and 7/3 (as well as Tq, and 7 q,) so that if T < min{T, 3 , Tq,} and 


It will prove useful to adapt the following alternative representation for the function "0: 

(3.4) lA(A) = ^R+{X)h{X), X G (a, (3), 

ZTTl 

where the function i?(A) is defined viai?(A)^ = (A—a)(A — /3), with i?(A) analytic inC\[Q;,/3], and normalized 
so that R{X) ~ A as A ^ 00 . The subscript ± in i?±(A) denotes the boundary value obtained from the upper 
(lower) half plane. 


Following ||], we define a function g{z) as follows 

r0 

(3.5) 


g(z) = / log(z - s) ■!/^(s)ds, z G C\(-oo,/3], 

^ Oi 


where for each s, log(z — s) is chosen to be branched along (— 00 , s), log(z — s) > 0 if z > s. 
We now list several important properties of the function g{z)'. 

• 1 g is analytic for z G C\(— oo,/3], with continuous boundary values g±{z) on (—oo,/3] . 

• 2 There is a constant £ such that for z G [a, /3], 

9++9- -V(z) = £, 

and for z G R\[q:, /3], 

g+ + g--V <e. 

Indeed, for z > (3, 

g++g- — V — l = 2g — V — l = — [ Rhds, 


(3.6) 

Jp 

and for z < a, 

g+-\-g- — V — £ = — f Rhds. 

J Oi 

• 3 For z G [a, /?], 

fd 

9+- 9- = ^TTI / V'(s) ds, 

J Z 

and this function possesses an analytic continuation to a neighborhood of (a, (3). 
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4 For z> (3, 
and for z < a, 


g+- 9- = 0, 


g+- g-= 27ru 

Using the function g{z), and the constant £, we define (cf. Theorem 5.57 in 0 ) M{z) via 



and define Mi as follows: 

• For z outside the upper and lower lenses, 


(3.8) 


Ml = M. 


• For z within the upper lens, 


(3.9) 


Ml =M 


( 1 

y—g-N{g+-g-) 



• For z in the lower lens, 


(3.10) 


Ml = M 


( 1 

ygN{g+-g-) 



(Recall that g+ — g- possesses an analytic continuation to such a lens-shaped region around {a, /3).) 

With these two transformations, we arrive at a matrix-valued function Mi, satisfying the following 
Riemann-Hilbert problem. 

The problem is to determine a 2 x 2 matrix valued function Mi satisfying 


( 1 .) Ml analytic in C\Emi- 
(2.) Mi=I+ C>(i), z ^ oo. 

(3.) Ml possesses continuous boundary values {Mi)± for z G Umj. 
(4.) (Jump relation) (Mi)+ = (Mi)_Vm, 

where 


(3.11) 

(3.12) 

(3.13) 

(3.14) 


Vm 


Vm 


Vm 


Vm 


I ^Nig++g.-V-i) 

0 1 

1 0 \ 

g-JV(g+-g_) J J 

1 0 \ 
gJV(g+-g-) I J ’ 


0 1 

-1 0 


for z G 


J , for z G S 4 U E 5 , 
for z G Si, 
for z G S 3 , 

S 2 , 
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3.3. The construction of a global approximation to Mi. The next step in the Riemann Hilbert ap- 

(A) (A) 

proach is to construct an explicit approximation to Mi , which we will denote M^ . The definition of M^ 
also is adapted directly from ||]-^. Let i?“ denote a disc of radius S centered at a, and let Bg denote a 
disc of radius 6 centered at /3. We will define outside Bg U Bg as follows: 


(3.15) 

(3.16) 


m[^\z) = 


7(z) = 


/ 7+7 

I 7-7 
\ -2i 

{z- (3)‘ 


7-7 

2i_ 

7+7 


(z — a) 4 

To define m["^\z) for z within Bg , we require the following auxiliary function 

2/3 


(3.17) 


j {-29{z) + y{z)+^f^^, 


z G Rf. 


As h{P) ^ 0, — {2g{z) — V{z)— £) behaves like (z — 0)^^^ for z near /3 (cf. relations (^)), and we may 
choose the branch so that $/ 3 (z:) is analytic in a neighborhood of /3, and 


(3.18) 


‘^ 0 {z) = c{z- f3)+0 (^(z - , c> 0, z^ p. 


Let the radius 6 of the disc Bg be chosen sufficiently small so that <i >;3 maps Bg injectively onto a 
neighborhood D' of 0. The decomposition of into four regions is as in ||^: = I U II U III U IV, as 

shown in Figure 3.19 below. 

C = ^/ 3 (^) 




/ \ I' \ 

II'\ 



my 0 

\ / TAP / 



D. 


c/3 


Figure 3.19. Decompositions of and D'. 


The four rays on the right of Figure p.l9| divide the C“Plane into four sectors, 0 < arg(C) < -f-, 


3 ^ 

arg(C) < 71, —TT < arg(^) < and < arg(()) < 0. The intersection of these sectors with D' defines 
the regions V, ifi, Ilb and IV^ (The precise angles separating the sectors are not important: we could, for 
example, replace ^ by any angle strictly between 0 and tt, etc.) Finally, the regions I, II, III and IV are 
defined to be the pre-images under d)(z) of I^, ib. III' and IV'. We define the contour "f^r in the C-plane 
to be the four rays shown in Figure p.l9 , with orientation inherited from the original contour Eg and the 
mapping 

We will need the following auxiliary functions in order to define the global approximation. 

Denote 


(3.20) 
and define 

(3.21) 


w := e 3 


:C\7,, 


C 


2x2 
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(3.22) 




Ai(C) Aiiu^C) 

Ai'iC) co^Ai'iio^) 
Ai(C) Ai(u;^C) 

Ai'iC) uj^Ai'iio\) 
—Lo'^AiiujC) 
Ai'iC) -Ai'iujC) 
AiiC) —uj'^AiicoC) 
Ai'iC) -Ai'iujC) 


Here (T 3 is the Pauli matrix 


1 0 
0 -1 


T -3 

, for C G /, 


, for C G //, 

1) 

, for C G III, 

b 

1 

, for C e IV. 


The next auxiliary quantity which we define is the analytic matrix valued function defined as 

follows. 


(3.23) 


Epiz) = 


7 

-iy” 


-7 

-zy 


$ 


0-3/4 


(At first blush, this matrix valued function is analytic on a domain which consists of a disc centered at z = /3, 
with the real axis deleted. However, a straightforward argument (as in ||^, Section 7]) shows that E extends 
to an analytic function on the disc Bg.) 

We are now ready to define within Bg : 

(3.24) m}^^(z) = E^iz) ^'^i^^iz)) for z G B^. 


We now proceed to the definition of the approximation m['^\z) within Bg . As the construction is 
entirely analogous to the construction of in 1^ , we will only present the formulae. More details on 


this construction can be fou nd in ||^, Section 7] , or |^, Subsection 4.4]. We define the transformation $a(z), 
analogous to defined in ( |3.17 ) above, 

/Q/y\ 

(3.25) i2gia)-Via)-£-i2giz)-Viz)-£))^^\, 

As hia) yf 0, 2gia) — Via)—£ — (2g(z) — Viz) — £) behaves like (z — a)^/^ for z near a (cf. relations (|3.6D), 
and we may choose the branch so that 4>q,(z) is analytic in a neighborhood of a, and 


z&B°g 


(3.26) 

^aiz) = —c(z — a) + O ^(z — a)^^ , c > 0, z ^ a. 
Next we define E^iz) as follows: 

(3.27) E^iz) = ) $-3/4. 

We now define the approximation within Bf: 

(3.28) m[^\z) = E^iz) («>a(z)) e2-34>r/3^3^ fo^. g Bf. 


3.4. The Riemann—Hilbert problem for the error. The next step in the Riemann-Hilbert approach is 
to define the ratio of Mi to m[‘^^ : 

(3.29) S = Miiz)[M[^\z)y\ zGC\Smi. 

It then follows that i?(z) satisfies the following Riemann-Hilbert problem: 

Problem 3.30. The problem is to determine a 2 x 2 matrix S'(z) satisfying 

(a) «S'(z) analytic in C\Ss, 

(b) «S'(z) possesses continuous boundary values 5'±(z) for z G Ss, 
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(c) S = I + 0{l/z), Z^OO, 

(d) S+ = S.Vs, zeSs, 


where Ss is defined in Figure 3.31. 



(3.32) 

(3.33) 

(3.34) 

(3.35) 

(3.36) 

(3.37) 


^ 5.1 = Ml 


(A) 


The jump matrix Vs is defined as follows: 

1 QN { g ^+ g -- V - 

0 1 

1 0 

g-iV(g+-g_) 


ys,2 = Ml 


(A) 


V(S.3 = Ml 


(A) 


r(A} 


zeEsT, 

ZGSs.2, 

gA(g!-g_) 

zeEs,4, 


I ^N(g++g.-V-£) 
0 1 


ys,4 = Ml 

Vs,5 = , z G Es, 


^"s .6 = , Z G Es. 6 . 

The fact of the matter is that 

1. On the upper lens, lower lens, and the real axis component of E 5 

(3.38) \\Vs-I\\l^, \\Vs-I\\L-<Ce-‘^^, C,d>Q. 

This can be deduced from (3.32)-(3.36), because (1) the global approximation m\^\ as well as its 
inverse, is uniformly bounded on these por tions of the contour, ( 2 ) the off-diagonal entries of the 
middle matrices on the right hand side of ( 3.32 )- (|3.36|) are uniformly exponentially decaying in N 
on these portions of the contour, and (3) for z ^ 00 , z G Eg^i U Es^, one has < 

Cg-dN{\z\+i) j-i-Qj. complete argument, see the proof of [||, Proposition 7.64]). 

2. On the circle Es^g, we have the following complete asymptotic expansion: 


(3.39) 

(3.40) 






5(z — a)^/^ 


72(z — /3)^/^ jp R{s)h{s)ds 
7(z-/3)1/2 


72(z — q:)i /2 R[s)h{s)ds 


(3.41) 


(3.42) 


35 

2 2592 


'/3 


R{s)h{s)ds 


-2 


— 1 i 
i 1 

1 i 

i —1 

-1 12i 

-12i -1 


More generally, we have for k even, k > 2: 


vr = (^ 


R{s)h{s)ds 


— k 


^fc) 

2 2 

-i{sk-tk) 

2 2 
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and for k odd, A: > 1: 


(3.43) 




2 ^= Hkiz-a) 




(z-/ 3 )i /2 R{s)h{s)ds 

^ 2'^-^Skiz-f3y/^ 

(z - a)i /2 R(s)h(s)dsJ 


— 1 i 


i —1 


where 


(3.44) 

(3.45) 


So = to = Ij 

r(3fc + i) 


— 


54^/c!r(A:+ |) 


6fc + 1 

TT’ = - - -Sk tor A: > 1, 


6/c — 1 


3. On the circle dB'^, we have the following complete asymptotic expansion: 


(3.46) 

(3.47) 


1 + Ar2 

5(z - / 3 )i /2 




= 


72(z — q;)i/2 R{s)h{s)ds 
7{z — a)^/^ 


+ 


72{z - / 3 )i /2 R{s)h{s)ds 
-2 


(3.48) 


= 


35 

2 ^ 


R{s)h{s)ds 


— 1 —i 

—* 1 

1 —z 

—i —1 

-1 -12i 

I2i -1 


More generally, we have for k even, k > 2, 

2k-l 


(3.49) 


= 


_ Sk + tk i{tk — Sk) 

Ris)h{s)ds)’" V Sfc + tk 


and for k odd, fc > 1 , we have 


(3.50) 


Vk^ = 


2'=-isfc(z-/3)i/2 


( 2 ;-a)i /2 R{s)h{s)ds) 
2'^-Hkiz - a)y^ 
(z-/3)i/2 R{s)h{s)dsy 


k \ —i 


— 1 —i 

1 


— 1 i 
i 1 


The reader may verify these asymptotic expansions using the definitions (|3.36| )-(3.37) of Vs, the definitions 
(3.24) and (3.28) of , together with the well-known asymptotic expansions for the Airy function and 
its derivative. In Q, the authors present a complete proof of this fact (see. Lemma 7.34 and its proof), and 
we leave the reader to carry out the analogous arguments in the present case. 


3.5. The solution S of the Riemann Hilbert problem 3.30 and its asymptotic e xpan sion. In this 
subsection we will present a formula for the solution S of the Riemann-Hilbert problem 3.3C , as well as an 
asymptotic expansion for S in powers of . To do so, we will require the integral operators C±, defined 
as follows: 


(3.51) 

(3.52) 


C± : L\Y.s) ^ L"(Ss), 


1 


/(s) 


C±{f){z)= lim -^tSiV-ds, ze^s- 

z'^z, zG ± Side of Es 27rz s — z 


The relevant integral operator Cvs ■ is defined via 

(3.53) Cvs{f) = C.[f{Vs-I)\. 
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The asymptotic expansions described in ( 3.39 )-( 3.48 ), together with the estimates ( 3.3S ), and the well 
known boundedness of the integral operators C± imply that the integral operator Cvs is bounded, with 
operator norm satisfying 


(3.54) 


ic. 


VsllL2(Eg)^L2(Sg) < 


c 

N' 


This implies that I — Cvs can be inverted by Neumann series for N sufficiently large. Therefore, we may 
define 


(3.55) 


fMs = {l-Cvsr^Cvs{I). 


Then, as in Theorem 7.70], the explicit formula for the solution S to the Riemann-Hilbert problem (3.30|) 


IS 


(3.56) 






'Ss 


ds, z G C \ Eg. 


Furthermore, as shown in Corollary 7.77 and Theorem 7.81 in |Q, S is uniformly bounded in C \ Es, and 
possesses a complete asymptotic expansion in powers of N~^ (in Section 7 of @], the error matrix is called 
R{z), whereas here the error matrix is called S). 

The terms in the asymptotic expansion may be obtained in a number of ways. For example, one may just 
compute the Neumann expansion of /rg, and keep only those terms which arise due to integrations around 
the circles dBg and dBg. Alternatively, one may posit an asymptotic expansion for S in powers of N~^, 


plug this expansion into the Riemann-Hilbert problem 3.3C, and obtain a sequence of Riemann-Hilbert 


problems for the coefficients in the expansion, which may be solved iteratively. The proof of Theorem 7.81 
in Q| involves (1) defining an integral operator associated only to the circles dBf and dB^, using only 
the first I terms in the expansion of Vs on these circles, (2) building an explicit approximation Si{z) by 

using the first I terms in a Neumann series for ~ (^) comparing this approximation to the 

true solution S{z). The ensuing Riemann-Hilbert problem for the error is shown to have a solution which is 
uniformly bounded by the first neglected term in the asymptotic expansion, etc. 

The result of these calculations is the following asymptotic expansion for the solution S: 


(3.57) 


Siz)r^I + Y,N-^S,i2 

1=1 


where Sj{z) is piecewise analytic, and can be defined iteratively via the following sequence of additive 
Riemann-Hilbert problems. 


Problem 3.58. The first problem is to determine a 2 x 2 matrix S'l)^) satisfying 

(a) 5'i(z) analytic in C\{dBg U dB^), 

(b) 5*1 ( 2 ) possesses continuous boundary values (5'i)j_ {z) for 2 ; £ dB"^ U dBg , 

(c) = 0 ( 1 / 2 :), z ^ 00 , 

(d) (5i)+ - iSi)_ =Vi, ZG dBf U dB^, 


Problem 3.59. Having determined Si, we now determine Sk, k > 2 satisfying 

(a) Skiz) analytic in C\{dBf U dBg), 

(b) Sk{z) possesses continuous boundary values (5'fc)j_ (z) for z G dBf U dBg , 

(c) Sk = 0{lfz), z ^ 00 , 

(d) iSk)+ - iSk)_ = Vk + EjZl iS,)_ Vk-J , z G dBf U dB^ 

Each of these Riemann-Hilbert problems possesses a unique solution, explicitly representable by contour 
integration. It is interesting and useful to note that in fact each of these contour integral representations can 
be evaluated explicitly. Thus, for example, we include an explicit representation for S'i(z) in Appendix A. 
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Some very important properties of the sequence {5'fe}fe>o are established in the following Lemma. We will 
need the Pauli matrices 


(3.60) 


= 


1 0 
0 -1 


= 


0 -i 
i 0 


(Ti = 


0 1 
1 0 


Lemma 3.2. For each k S Z>o, there are complex valued functions Sk{z) so that the following holds. 

1. For k even, the function Sk{z), defined to be the unique solution to the Riemann-Hilbert problem 
3.5i satisfies the symmetry conditions 


(3.61) Sk{z) = s^^\z)I + s)^’{z)a2. 

2. For k odd, the function Sk satisfies the symmetry conditions 


,( 2 )/ 


(3.62) Skiz) = sl^\z)a3 + slf'(z)ai. 

The functions sl^\z), j = 1,2, are piecewise analytic functions of z, explicitly computable in terms of the 
parameters t. 

In words: for k odd, Sk{z), is a symmetric, trace 0 matrix. For k even, Sk is the sum of a multiple of the 
identity matrix, and a skew-symmetric matrix. 

Proof. For fc = 1, the Lemma is seen to be true by studying the solution to the Riemann-Hilbert problem 


o(2)/ 


3.58 


(3.63) 




-ds. 


s — z 


where the i nteg ral is t aken along the boundaries of the discs i?“ and Bg, oriented clockwise. Since Vi, as 
defined in (3.40) and (3.47) is clearly symmetric trace zero, inherits the same property. 

We now argue inductively. Suppose that k is even, and the Lemma holds true up to fc — 1. The reader 
may verify that the right hand side of the jump condition (d) in the Riemann-Hilbert problem 3.59 is of the 
form 


(3.64) 


ui^\z)I + u]f>{z)a2 


( 2 ), 


where v^^\z), j = 1,2, are scalar complex valued functions of z. Indeed, for k even, Vk{z), as defined in 


(3.42) and ( 3.49D are of this form. Moreover, each term appearing in the sum on the right hand side of the 
jump condition (d) in the Riemann-Hilbert problem p.5£ is also of the same form (this can be deduced from 
the form of Sj and Vk-j for each value ofj € {1,...,A: — 1}, given that k is even). There is an explicit 
formula for the solution to this additive Riemann-Hilbert problem, which is gi ven b y the Cauchy transform 
of the right hand side of ju mp c ondition (d). Since every term is of the form ( 3.64 ), the Cauchy transform 
is also of this form, and so ( tl. 6 l| ) holds true for k even. 

The argument for k odd follows the same line of reasoning, and we leave the details to the reader. 

□ 

At this point we have established the following formula for the solution Y to the original Riemann-Hilbert 
problem stated just before Subsection 2.1: 


(3.65) 




where I is the Lagrange multiplier associated to the variational problem, Slz) is the solution to the Riemann- 
Hilbert problem 3.30 for the error, M^^\z) is the explicit approximation defined in ( 3.15 ), ( 3.24 ), and ( ^.28 ), 
T{z) is either the identity matrix I for z outside the lens shaped regions, or one of two explicit triangular 
factors appearing in (3.9) and ( |3.10 ), and g{z) is the function defined in (3.5). It is straightforward to verify 
that with the exception of S(z), all terms appearing in ( 3.65| ) depend analytically on the times t for t in 
any open subset of T(T, 7 ), and in fact extend as infinitely differentiable functions to all of T(T, 7 ). It is 
also true that the function S(z) enjoys the same analyticity properties. This is made precise in the following 
Theorem. 
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Theorem 3.3. The matrix valued function S is infinitely differentiable in the times t for t G T(r, 7), 
and derivatives of S again possess an asymptotic expansion in N, obtained by differentiating the original 
asymptotic expansion for S. The coefficients in this asymptotic expansion depend analytically on t for 
t G T(T, 7) (with smooth extension to the closure), for T sufficiently small. 

Proof. The proof hinges on the following fundamental facts: the jump matrix Vs is analytic in a vicinity of 
every point z on the contour Ss, it possesses an asymptotic expansion in powers of N~^ in which each term 
is a piecewise analytic function of z and depends analytically on t for t G T(T, 7); this asymptotic expansion 
may be differentiated term by term with respect to the times t G T(r, 7), yielding asymptotic expans ions for 
arbitrary derivatives of Vs- From this information, it is straightforward to prove that p,s defined in ( t 3 . 55 | ) is 
an analytic function of t for t G T(r, 7), has an analytic continuation to a neighborhood of the contour Eg in 
the z-variable, and possesses an asymptotic expansion (which can be differentiated term by term) obtained 
directly from a Neumann series representation. An outline of the argument is as follows. 

First, observe that the terms in the Neumann series representation of p,s ^re individually analytic functions 
of z in a vicinity of every point of the contour E5, and they are analytic function of t for t G T(T, 7). Second, 
the series is uniformly convergent. This follows because ( 1 ) each Neumann iterate can be expressed as an 
integral on a deformed contour which is bounded away from the contour Eg and (2) the j-th iterate may 
then be bounded by C{D/Ny, where the constants C and D are independent of j. For derivatives of /is, 
similar calculations, using the fact that derivatives of Vs possess asymptotic expansions, show that an Lth 
order derivative of the jth iterate may be bounded by C{D/Nfi, where C and D are independent of j, 
but may depend on L. This shows that the Neumann series for /is may be differentiated term by term. 

The proof that derivatives of S possess expansions which may be obtained by differentiating the original 
expansion for S term by term now follows from the uniform convergence of the Neumann series for /ts and 
its derivatives. □ 


4. The mean density of eigenvalues: exact formula and asymptotics 


The previous section explains (and extends to higher orders) the work in |^-Q, in which techniques for the 
asymptotic analysis of Riemann-Hilbert problems, originally developed for the analysis of singular limits of 
integrable partial differential equations by Deift and Zhou (and further developed by many researchers), are 
used to establish an asymptotic description for the polynomials pn{z) and pN-i{z) orthogonal with respect 
to the A^-dependent measure ( | 1 . 23 | ). There are many consequences of the complete asymptotic description for 
the orthogonal polynomials. For a number of such applications, we refer the reader to where detailed 

asymptotic expansions for (1) the leading coefficients of of the polynomials, (2) the recurrence coefficients, 
( 3 ) the zeros of the polynomials were obtained, along with a proof of universality of spacings of eigenvalues of 
random matrices for a wide class of unitarily invariant probability measures on Hermitian matrices. In this 
section our goal is to establish formula (^), and prove that this formula for the mean density of eigenvalues 
holds true for z in a fixed si zed neighborhood of [z*,(f\ for N sufficiently large. The parameter z* will be 
chosen below. The f ormu la ( 4 N) holds true in a large neighborhood containing the point / 3 , and there is an 
analogous formula, ( 4 . 15 ) which holds true in a fixed size neighborhood of the interval [a,z*], and the two 
formulae hold simultaneously on an overlap region. In order to achieve this overlap region, we will need to 


introduce T,a, 7/3, Ta, and 70,, which are values of T and 7 from Theorem ^ so that if T < imn{Ti3,Ta} 


and 7 > max{7^,7a}, then the neighborhoods I?e,/3 and may be taken large enough to have nontrivial 
overlap. _ 

Using formula ( | 4 . 4 D we will obtain a complete asymptotic expansion for the mean density of eigenvalues, 
and then in the following section we will evaluate an asymptotic expansion of integrals of the form 


( 4 . 1 ) 


f{\)p%\X)dX, 


for general C°° functions that grow no worse than polynomial at 00. 
We begin with two important observations. 
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Observation 1: It is possible to use the results of the previous section to compute an asymptotic 
expansion for the mean density which is valid for all A G {a,P). That expansion is as follows, 


(4.2) 


p«(A)= m + ^ ■=»{"/ '/'(»)* 


+ 


Ar2 


H{X) + G(A) sin ^Pis)ds 


in which H (A) and G(A) are locally analytic functions which are explicitly co mput able in terms of the original 
external held I^(A). However, for the computation of integrals of the form (^), this expansion provides a 
piece of the puzzle, but clearly cannot be used in a vicinity of the endpoints a and/3. Indeed, the asymptotic 
expansion fails in a vicinity of these endpoints, because the second term in (O) possesses poles, whereas 
the mean density does not. From this observation we see that to compute asymptotic expansions of 
integrals of the form ( [4.1[ ), we will need to to have an asymptotic expansion for which is valid for A near 
these endpoints. It is a piece of good fortune that the representations which we obtain are actually valid on 
regions which cover the entire interval [a,/3], and we will not need to evaluate integrals “in the bulk”. 

Observation 2: It is straightforward to use the results of the preceding section to prove that p^^(A) is 
exponentially small (in N), and exponentially decaying (in A) for A G R. \ (a — <5, /3 + J), for any fixed 5 > 0. 
We will not present the details of those calculations, but leave them to the avid reader. The implication of 
this is that the contribution to an integral of the form (4.1) from the set R.\ {a —5, (3+S) is exponentially small 
in N, and hence negligible to all orders in N. Therefore we must only compute an asymptotic expansion for 
integrals of the form 

(4.3) / /(A)pW(A)dA, 


/ oc — S 


in which the allowable functions / are C°° smooth, and compactly supported within (a — 6, (3 + S). 

We now turn to the formula (^^) for pj^^, and re-express this fundamental quantity using the Riemann- 
Hilbert analysis of the previous section. We present two functions, and representing pj^^; the first 

one is a valid representation for p^^ in an interval containing z = /3, and the other, in an interval containing 
z = a. The first such formula is the following explicit, exact expression: 

(i-i) = (^y - |2 Ai (<!.«(=)) (4.M.))| 


( 24 *'(d>p( 2 :))) -^piz) (^Ai{^p{z)) 

(^'H4-(4>p(z)))^^ (5H 4-(4>p(z)))^^ 

- (5H vl-(4>p(z)))^^ ( 5 ' H vl/($p(z)))^^' 


where 


(4.5) 


vh(C) = 


r { Azio 

A^MC) \ 

J V At'iO 

u^At'MO J 

] ( MO 

—uOAi{LoC,) \ 

[ V ^*'(0 

-Ai'{ujQ ) 


for ( G C+, 
for C G C_. 


The transformation $p(z) was defined in (3.17), 7 ( 2 ) was defined in (^.16|), and the matrix B is defined by 


(4.6) 


B{z) = 


7 H^) -7(2) 

-Z 7 “^(z) —ij{z) 




0-3/4 


and S{z), defined in ( 3.29| ), is the solution to the Riemann-Hilbert problem ^.30| for the error, which possesses 
a complete asymptotic expansion (3.57), with coefficients Sk solving the Riemann-Hilbert problems 3.5S and 
3.59, and satisfying the symmetry conditions described in Lemma 
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Shortly we will present an analogous representation of near the left endpoint a. The following lemma 
will be instrumental in describing the relation between these two representations as well as their respective 
domains of validity. 


Lemma 4.1. There is a unique point z* in the interval (a,/3) at which ($/ 3 ( 2 :*)) = ( 4 ) 0 ( 2 :*)) 
Proof. It follows from formulas (L17), ( 3.25D and ( |3.6| ) that for z < (3, 


(4.7) 

and for z > a, 

(4.8) 




3N 




-hds; 


= (™) I 


It then follows from comparison to ( 3^) that the integrals in ( |4.8[ ) and (4.7) are both positive. 

By construction (see ( 3.18 ) and ( 3.26 )) we know that ^is{z) and ^a{z) are negative real for a < z < (3. 
The z* we seek must satisfy 


Rhds ± 


Rhds = 0. 


We already know that 


Rhds - 


Rhds = / Rhds = 27ri. 


It follows that z* is uniquely determined, equivalently, by either of the two equations: 

rd 

Rhds = 771 = 


J 2* 


Rhds 


□ 


Lemma 4.2. There is a Tp > 0 sufficiently small and 7 ^ sufficiently large (both independent of N) so that 
fort G T (Tfl , 7 / 3 ), there exists a complex neighborhood B^g of the interval [z*,2], independent of N, so that 
formula (i4-4 ) holds true for all z G Up. 

Proof. The proof is a straightforward, though lengthy algebraic manipulation of the explicit relationships 
between the original quantity Y(z), and the error matrix S{z). The only analytical issue is this: is the 
transformation ^p{z) an invertible analytic transformation from some fixed size neighborhood of the interval 
[z*,2] into its range, for t sufficiently small? But this is clear, because for t = 0, the function 4)g is quite 
simple, 

2/3 


( 3JY _ _ \' 

— Vs-2Vs + 2 dsj 


The reader may easily verify that there exists a branch of this function which is analytic on a fixed open 
neighborhood of the interval (—2,2]. Now since ^p(z) depends analytically on t, it is clear that ^p{z) is 
certainly an invertible analytic transformation from some fixed size neighborhood of [ 0 , 2 ] into its range. 

The reader may verify that formula (4.4) is true, by starting with the formula (O) for and re¬ 

expressing this formula in terms of the solution Y of the original Riemann-Hilbert problem 

(4.9) p[^'^ {z) = Ki ( 2)^21 [z) - Tn {z)Yi, ( 2 )] . 

Next, observe that for 2 G C± n {lens shaped region in Figure 0 } , we have the following representation 
for the solution Y of the original Riemann-Hilbert problem: 

Y{z) = 

e^^^S{z) Ep{z) ( ±eT^L-.-) J 

= e^'*35(2) Ep{z) 4'($/3(2)) )/6)-3^ 


^N{g- 


(4.10) 
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Here Ep{z) was defined in ( 3.23| ), and is related to the matrix valued function B{z) defined in (4^) via 
(4.11) Ep{z) = ^E^/^B{z). 

The simplification in (4.10|) follows from the identity 


(4.12) 


1 0 
Tl 1 


^ ) „N{g-^)r73 _ §V{z)a3 

±eT^(9+-9-) 1 ' “ 


which can be deduced from the properties of g(z) and ^/ 3 {z). Next, from the definition (O), one easily 
obtains the identity 


(4.13) 




The following identity for 4' is also equally easy to prove, using the definition (^^) and the differential 
equation satisfied by the Airy function: 


(4.14) 


' I ''(0 = ( ^ J ) '*^( 0 - 


Formula (4.4) may now be verified by differentiating (4.1C), substituting into (4.9) and using the differential 
relations (4.13) and (|4.14). 

□ 


Formula (4.4) may seem unwieldy due to its length. However, this is not quite the case. Indeed, the 
first two lines of this formula do not contain S(z) at all, and are explicit combinations of Airy functions 
(composed with ^/ 3 {z)) together with other functions defined explicitly in Section ||. Moreover, we will show 
that the third and fourth line of (4.4) possess asymptotic expansions in powers of N~^ in which each term 
is of the same explicit form as the first two lines. This is the content of Lemma 4.4 below. 


There is, of course, an analogous formula which is valid on a fixed size complex neighborhood of the 
interval [—2, z*]. We summarize the analogous result in the following Lemma. 

Lemma 4.3. There is Ta > 0 sufficiently small and ■ja sufficiently large (both independent of N) so that 
fort G 'T(Tq,, 7 q,), ther e exists a complex neighborhood of Ha of the interval [—2,z*], independent of N, so 
that the formula (j.lL) holds true for all z € Bq,. 


(4.15) 


f , l'{z) 


W“^(z) = - 


\4$a(z) 7(z) 


[2 Ai{‘^a{z)) Ai'{^a{z))] 


where in this case 
(4.16) 


-<^'a{z) (Ai'($«( 2 :))) - ‘t>a{z) (^Ai{<t>a{z)) 

+ ^ [(5' B 4/ (4>„(z)) ) (5 H 4/ (4>„(z)) ) 

- H ^ (<^a(z)) ) ( 5 ' B 4- ($„(z))) 

B(z) = ( 

ij(z) 17 "(z) J 


Next, we will compute an asymptotic expansion for by using formulae (^), ( 4.15| ), and the asymp¬ 
totic expansion for S(z) described in Section^. 
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Lemma 4.4. If t G T(r^, 7 ^), then there is a fixed size neighborhood of the interval [z*,0\, the following 
asymptotic expansion holds true: 


(4.17) 


I 

m L 


( 5 ' B B ^($^(^)) 

- (^S B 4'(4>/3(2)) 


= E 

j even, j>2 


(s' 


j even, j >2 


7(z)2 

nA^pjzfMzf 


' E 

j odd, j>3 


N-^dfiz)'^ii{^p{z))^2i{'I’0iz)), 


where for each j, the coefficient functions aj(z), bj{z), and cfiz) are jointly analytic in z and t, for z in a 
fixed sized complex neighborhood of the interval [z*, /3], and for t G ir(T/ 3 , 7 , 3 ). 


Proof. From Lemma 3.2, we have an asymptotic expansion for the function S{z): 

C/„\— T , \ ' I I AT-k 


S{z)=I+ Y. 

k odd, k>l 


+ 


E 


k even, k >2 


(^s[^\z)I + S^^\z)a2^ 


N 


— k 


in which (z) are functions of z that are piecewise analytic functions of z, which are, in particular, 
analytic on a fixed size neighborhood of the interval [0,/3], and de pend analytically on t for t G T(T 3 , 7 , 3 ). 
An asymptotic expa nsio n for S'{z) is obtained by differentiating ( [l.l8[ ) term by term. This, together with 
the explicit formula (4T) for B{z), may be substituted into the left hand side of ( 4.17 ). Now straightfo rwar d 
algebraic manipulation yields an asymptotic expansion of the form displayed in the right hand side of ( 4.17 ). 
The quantities dj and bj may be computed via these calculations, but all that we will require is that these 
functions are analytic in a neighborhood of [0,/3], and depend analytically on t for t G TiTp,"//}), which is 
manifestly true. □ 


There is, of course, an analogous result which holds for the interval [a — 6, z*], which we will only state, 
as the proof follows by the same arguments as the proof of Lemma 4.4 

Lemma 4.5. //t G T(Tq, 7 q), then there is a fixed size neighborhood of the interval [a,z*], on which the 
following asymptotic expansion holds true: 


(4.19) 


% 

Yn i 


( 5 ' B ^ ($„(z)) ) B T ($a(2)) 

-( 554 - ($„(z)) ) ( 5 ' B T ($„(z)) ) 


2 IJ 


J even, j>2 

+ Y iv-^'6,(z)^i%ML 

.even, ,>2 

00 

+ Y Af-^gj(z)4'ii(4>a(z))T2i($o(z)), 
j odd, j>3 


where for each j, the coefficient functions d. (z), bj{z), and cfiz) are jointly analytic in z and t, for z in a 
fixed size complex neighborhood of the interval [a,z*], and for t G T(Tq,, 7 q). 
















30 


N. M. ERCOLANI AND K. D. T-R MCLAUGHLIN 


Remark: Observe that if we take T < min{r^, Ta} and 7 > max{ 7 ^, 7 ^} then the conclusions of Lemmas 
4.2, 4.3, 4.4, and 4.f hold true simultaneously. From now on we will take t S T(T, 7 ) for such values of T 
and 7 . 


5. Proof of the Main Theorem 


We are now in a position to analyze the structure of the fundamental moments ( 1.33| ) which will enable us 
to establish the 1/iV^ expansion of logZff. To do this we will be using the representations of the one-point 
function developed in the previous section. Precisely, we use 

(5.1) Pn\z) = XaWPA + Xp{z)Pn'^'^ 

where {Xa,Xp} is a partition of unity for K defined by 


(5.2) 


Xpiz) is 
0 < Xpiz) < 
suppxp C 
Xp{z) = 1 


C°° with 

1 

{z* - e,oo) 
for z G {z* 


We then define Xa by 


Xa{z) = 1 - Xp{z) 


which has similar properties to xp but with support in (— 00 , z* + e) 

The expression (P) is merely a different representation of the one-point function , and the partition of 
unity allows us to carry out our asymptotic analysis separately on Xaiz)p^^’‘^\z) in Bq, and on xpiz)P^ n^\z) 
in Bp, and then, at the end, to add the two contributions together. This partition of unity will prove especially 
useful in avoiding the evaluation of boundary terms during repeated integrations by parts. Because of the 
complete similarity of the analysis for each of these separate pieces it will suffice to demonstrate our claims 
for xp{z)p^j!f’'^\z). 


We now begin the analysis of integrals of the form 
(5.3) 


0-\-S pz*-\-e nf3-\-S 

/(A)pW(A)dA= / Xa(A)pi^^(A)dA+ / xp(A)pW(A)dA. 

i — S Jot —6 Jz* — € 


(Recall from Observation 2 and formula (O) that we need only consider integrals over the bounded set 
{a — 6, P + S) for some sufficiently small (but independent of iV) (5 > 0.) We will prove that if the function 
/ is C°° smooth and compactly supported within {a — 5, P + 5), the second integral in (^.3|) above possesses 
an asymptotic expansion in even powers of N: 

^/3+(5 

(5.4) / xp(A)pW (A)dA = /o + iV-Vi + + ■■■ . 

Jz*-e 

As we have explained, this will complete the proof that 

(5.5) 


A^logZjv(t) =eo(t) + ^gi(t) + . 


i.e. that dti logZjv possesses a complete asymptotic expansion in even powers of N. Since this asymptotic 
expansion is uniformly valid for all t G T(T, 7), we may integrate (5.5) from t = 0, and thus we will complete 
the proof of Theorem 0- 


From the analysis of the previous section, in particular lemma 4.4, one observes that the terms appearing 
in an expression of the form X^xpWPn'^'* W ^be following four types of integrals: 
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1 /*/5+<5 

(O') J ^ 9{z)Foi'^0Wfidz, 

/•/3 + (5 

( 1 ') N~^ / g{z)Ai{^i 3 )Ai {^ 0 )dz {j odd), 


€ 

. /$« /•/ 3+<5 

) / g(2;)Az2($/3)(iz ( j even), 

Jz*—e Jz*—e 

ry{ 2 S^ , . 

/ / Kz)WY{^0)dz, (j even), 

Jz*-e ^/3’/3 Jz*-e 


(20 

(30 N-^ 


where g{z) and g{z) are general infinitely differentiable functions of z, and are compactly supported within 
{z* — e,P + S). The function Fq appearing in (OO above is defined by 


(5.6) 


Fo(C) := [AifiO-CAiHC). 


The main result would be established if we could show that all of these quantities possess asymptotic 
expansions in even powers of N. We therefore must analyze the following four quantities: 


(0) 

^ p(3-\-6 

N J . 9 {z)Foi<^0)^'pdz, 

(1) 

|*/3+(5 

/ g{z)Ai{^is)Ai' {^0)dz, 

J z* —e 

(2) 

rp+S 

iVi/3 / ~g{z)Ai^i<^>0)dz, 

J z* —e 

(3) 

r0+S 

iV-1/3 / g{z)[Atfi<^>0)d^ 

J 2* —€ 


and prove the following statements true: 

(0), (2), and (3) have asymptotic expansions in even powers of N. 
(1) has an asymptotic expansion in odd powers of N. 


In fact, once we show that integrals of the form (0) have asymptotic expansions in even powers of N, the 
following three basic observations take care of the remaining types of integrals. 

Observation 1: 

(1) has an asymptotic expansion in odd powers of N if and only if (2) has an asymptotic expansion in 
even powers of N: 


2 



gAi{^ 0 )Ai'{^f})dz 




Ai^{^l3)dz 



Ai^(^l})dz 


where in the first equality we’ve integrated by parts. (No boundary terms are present as the integrand 
is compactly supported within {z* — e,/3 + d).) We now observe that is analytic and nonvanishing 

throughout the region of integration so that the integral is well-defined. Lastly, the notation (j)'i^{z) := 


N ^/^$o(z) defines an analytic, nonvanishing function which is independent of N. The quantity \g/4''ii 


is 


a new infinitely differentiable function of z, compactly supported within (z* — e, f3 + 6 ), and we have shown 
that integrals of the form (2) possess asymptotic expansions in even powers of N if and only if integrals of 
the form (1) possess expansions in odd powers of N. 

Observation 2: 

Integrals of the form (2) have asymptotic expansions in even powers of N if and only if integrals of the 
form (0) have asymptotic expansions in even powers of N: 
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1 

rd+S 

IN 

N , 


1 

1 _ 


/ gAi^{^p)dz = 

Jz*-e Jz*-€ ^8 

where once again we have integrated by parts in the second equality. Now we have arrived at an integral of 
the form (0), since the quantity is infinitely differentiable on [z* — e,/3 + i5]. 

Observation 3: 

Finally, if integrals of the forms (0) and (1) respectively have even and odd asymptotic expansions, then 
integrals of the form (3) have even asymptotic expansions: 


rd+S 

/ g[Aif{<i> 0 )dz = 

J z* —e 


-N- 


t-d+s 


U-e 

1 

1 _ 


{\Ai{<^d)^i\^d) + ^^dPo{^d)}dz 


= -3^- 




rd+S 

/ 2 * _e 

i-d+s 


<P'd 


^'d 


Ai{^ 0 )Ai'{^ 0 )dz 


'^F 8 ^ 0 )^' 0 dz. 

90 


There is an integration by parts in the first equation. The two terms in the last equation are integrals of the 
form ( 1 ) and ( 0 ), respectively, and this establishes our third observation. 

We turn now to the heart of the matter. We will give an inductive proof that 

r/3+<5 


(5.7) 


1 fPA -0 

N y. 9 {z)Foi^ 0 )<^> 0 dz 


has an asymptotic expansion in even powers of N. We will integrate by parts repeatedly, peeling off con¬ 
tributing factors as we go along. The proof proceeds in four steps. The first three steps involve setting up 
an integration by parts inductive argument. 

Step 1: Show that the function Fq, which is analytic in C, possesses the following asymptotic expansion 
for C ^ —oo: 

(5.8) FoiO = (co + ci(-C)-^ + C 2 (-C)-'‘ • • •) 

(c?o(—C) ^ + dii—C) ®+d2(—C) -) sin ^-(— 

+ (/o(-C)-' + /l(-C)-" + /2(-C)-" + • • • ) cos (^(-C)'/") • 

Step 2: Show that Step 1 implies that there are functions Gi(C)) G 2 (C), and G 3 (C) which are analytic in 


C, satisfy 


(5.9) 

G'i(C) 

(5.10) 

G^(C) 

(5.11) 

G'siO 
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and possesses the following asymptotic expansion for ^ > — oo: 

(5J2) G,(0 = + + 

+ (/i‘'(- 0 -= + + - ■) ™ ( 5 (-C)’''“) 

(5^13) G,(0 = 

+ + ■-■)»» (5(-0“'“ 

+ + /P>(- 0 -= + - ■) ™ 

(5.14) G3(0 = *(-<)““''“ + ■■■ 

+ {C(-o-'+ /f >(-o-' + - ■) ™ (|(-o’/=). 

This defines antiderivatives of Fq which we will also need in what follows: 


(5.15) 

Fi 

- r 1 

- Gi -—1 

(5.16) 

F2 

- G 2 + —1 

(5.17) 

F3 

= G3-^ 

^ 105 


and we will show that these three anti-derivatives are all exponentially decaying for C ^ -foo: 
(5.18) 

Step 3: We express the integral (|5.7D as follows 




(5.19) 

(5.20) 

(5.21) 


^ ^/3+(5 

— / g{z)Fo{<^>0)^'f^dz = eo +A, 

iV 7^*—e 

" ‘^'pdz 

J ^ 5 ( 2 ) [Foi^p) - Co ^ j •^'pdz +j^J 9{z)FQ{'^0)^'pdz 


Co = 


r0+s 


Now Co is in fact independent of N: 

(5.22) eo = diz)-^ 0 (-^/3(2:))^'^^^ dz = —^ g{z)R+{z)h{z)dz 


and so we must show that A possesses an asymptotic expansion in even powers of N. Integrating both terms 
in ( 5.21 ) by parts, we obtain 


1 rP 1 

(5.23) A =J ^ g'{z)Gi{^p{z))dz - — g'{z)Fi {^^{z)) dz. 

Observe that in (|5.23| ), the boundary terms from z* — e and /3 -I- 15 do not contribute because g is compactly 
supported within (z* — e, P + 5). Furthermore, the boundary terms from the integrals &t z = j3 cancel, 
because Fi($^(z)) and Gi($^(z)) coincide a,t z = (3. Since Gi and Fi are both uniformly bounded, it 
follows immediately that < GN~^. However, it is not yet straightforward to deduce the asymptotic 
behavior for H as ^ 00 , and that is our next task. 
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Two more integrations by parts yield 

1 rf^ 

( 5 ^ 24 ) 41 = -- 




G^(‘if{z))dz 


1 

'iv 


/3+5 






F 3 i<^f3{z))dz. 


Recalling the asymptotic expansion (5.14) for G 3 valid for C. — 00 , we may rewrite (5.24) as follows: 
(5.25) A = eiN-'^ + A 2 , 


(5.26) 

(5.27) 


ei = -N 


Ao — — — 


N 


k [kH 
k [k^i 


k 


r0+s 


k [k/^k 


I 

{G3{^i3{z))-cf (-$^(z))l/2) 
F 3 i^niz)) dz. 


As with Co, Cl is independent of N. Indeed, using (|4.7D, we have 


(5.28) 

(5.29) 


ei = - 


i: 


k {kA\ 


Co i-(j)f3iz)k^‘^ dz, 


where for z G {z* — kpiz) = — 


-3z 


2/3 


R+{s)h{s)ds 


Observe that once we have established Steps 1 and 2, then Step 3, carried out above, shows that 


(5.30) 


1 1 

— J ^ g{z)Fo{^ 0 )(^>'f^dz = eo + ^ei + A 2 , 


with Co and Ci independent of N (see (5.22) and ( p.2^ ), respectively), and with A 2 defined in (5.27). 
Furthermore, it is straightforward to prove that IA 2 I < and hence ( 5.30| ) shows that the first two 

terms in the asymptotic expansion of (5.7) are in even powers of N, as desired. The above calculations must 
now be automated to establish the existence of a complete asymptotic expansion in even powers of N. 

Step 4- The procedure should now be clear to establish the induction. We must prove that we can 
integrate by parts thrice, peel off the next term in the asymptotic expansion, and repeat. 

We must prove that there exists Gm, Fm, m = 4, 5,..., satisfying a number of relations, all appearing in 
triples. First the differential relations: 


(5.31) 


Gk+i = G 3 j - c^’(-C)'/^ G' = G3J+1, G' = G 


43i+2, 


(5-32) F^+i - F 3 j, A3j_|_ 2 - Tsj+i, ^( 3^+3 - ^ 3 ^+ 2 , 

then the asymptotic prescription 

(5.33) G 3 J +3 = (-C)'/"(4''+'^+cl"^+")(-C)-' + ---) 

4in (5(-C)='4) + Gr«' cos (|(-C)»'4), 
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where and are asymptotic series which behave as follows: 

(5.34) 


(5.35) 


Gf+3) ^ (_^)-(3,+5)/2 j'^(3,+3) ^ ^(3j+3) ^ ^ 

For J even: { ^ ^ ^p,+ 3)^_^^_3 ^ ^ ^ ^ ^ 

(^^3i+3) ^ ^_^^-(3 j+8)/2 |'^(3j+ 3) ^(3j+3) ^ 

For j odd: { ^{ 3 j+ 3 ) ^ ^_^^-( 3 j+ 5)/2 ^ ji 3 j+ 3 ) j{ 3 j+ 3 ) (_^)-3 . . . ^ 

(note that this implies a similar asymptotic behavior for Gs^+i and G 3 j+ 2 ). The quantities Fg must satisfy 
the following estimates for C ^ + 00 : 

(5-36) iTsj+il, |T3j+2|, lAj+sl < G3je~3‘^ , 

and finally we require the following relation between and Gf. 

(5.37) G™(0) = F„(0), m = 3j + 1, 3j + 2, 3j + 3. 


Once we have the existence of the functions Fi and Gi out of the way, then we may proceed with the 
inductive integration by parts. Suppose we have arrived at the situation 

2 rP+S 

(5.38) — J g{z)Fo{^0)^'i3dz = io + N '^ei + --- + N 


where Cj are all independent of N, and 

(5.39) /(.-)) 

^(G3i{<^0{z))-cf^\-^0Y/^) dz 


(_ 1 )£ rf3+5^ 


N 


Ci (w‘I: ((-y-■ 4 oy 

>^F3l{^l3)dz, 


and in each nested set of derivatives appearing in the integrals in ( ^.39 ), the differential operator 

iMF) 


(5.40) 


appears 3£ — 1 times . We further assume that in ( |5.39D , G^i satisfies ( 5.33| ), F 31 satisfies ( |5.36|) , (bot h wit h 
j = £ — 1) and (5.37) holds with m = 3f. In order to integrate by parts, we note that from (5.31) and (^.32), 


(5.41) 

(5.42) 

Using 

(5.43) 


{G3,{^p{z))-cf\-^pY/^) = ±Ag3,+i($^(z)) 


F3i{^l3{z)) — ,, (^x —F3t.+l{^l3{z))- 


^'^{z) dz 

and integrating by parts, we find 

^r+i = 

(_1)^ rd+s d 


N J^._^dz 
(- 1 )^ 


N 


,dz 






X {Gse+i) dz 

xF3e+i{^l3)dz, 
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where now the differential operator {d/dz{^'p) ^ •) appears times. At this point, using (i3)> it is straight¬ 
forward to prove that |Af+i| < . However, we may integrate by parts two more times, and we 

then arrive at 


(5.44) 


A^+i = 



X (Gse+s) dz 

■■A((5.irsu))-))x 

xF3t+3{<^0)dz, 


where now the differential operator ^ 

may now be rewritten as follows: 


•) appears 3£ + 2 times in each integral. Expression (5.44) 


(5.45) 

(5.46) 


(5.47) 


Af+i — N ^e^+i + A^+2) 

£i+i = 

O(‘■"A- s(w-‘ X 

X (4”"“' (-tri''") d.. 


Ai+2 = 



xF3e+3{<Pl3)dz, 


where the differential operator {d/dz{^'p) ^ ■) appears 3£ + 2 times in each integral. The quantity e^+i is 
easily seen to be independent of N using the representation (g3) of $, 3 : 


(5.48) 


e^+i = 



X dz. 


Therefore we have deduced that 
(5.49) 


j-lB+S 

— I g{z)Fo{^f3)^'f3dz = eo + N~‘^ei -\ - \- + A 1 + 2 , 


with e^+i defined in ( 5.48| ), and A ^_|_2 defined in (5.47). Summarizing, once we establish the existence of 
the functions Fg and Gi, we have established by induction that (5.7) has an asymptotic expansion in even 
powers of N. 

Proof of Step 1: We establish Step 1 using well-known [|p asymptotic representations of the Airy 
functions for large negative values of the argument: 
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(5.50) 

Ai{—z) = 

1 

0Fzi/4 

(5.51) 

Ai'{—z) = 

zl/4 / 

V 

where 



(5.52) 




(sin(C + ^)fi'i(C) - cos(C + ^)W 2 (C)) , |arg 


, 2 
z|<-. 


|arg z| < -77, 


C = 

and the Wi{Q are asymptotic series in large C defined recursively in terms of the coefficients 

(2s + l)(2s + 3)(2s + 5)... ( 6 s - 1) 


Us = 

6 s + 1 

Vs = -T- tU. 


(216)«s! 


with 


6 s- 1 

Mo = 1, "^0 = 1- 


Explicitly the asymptotic series are given by 


(5.53) 

s^O ^ 

W2 = 


s^O 

(5.54) 

oo 

T\sV2s 

w3 = 2^{-i) 

s=0 ^ 

oo 

W4 = 

s=0 


, 7'-2s+1 
^ 2 s+l ’ 


We caut ion the reader that in formulae (|5.50D -( 5.54|) the v ariab le z > 0 is a dummy variable, as is C appearing 
in (|5.52| ) above. We will replac e z ap pearing in (|5.50| )-( 5.54 ) with —^ (where C, < 0). Consequently, the 
parameter Q appearing in ( 5.50 becomes (2/3)(— 


One can now substitute the expansions of the Airy functions (p.50|) into the expression (^) for Fq{Q 


and collect separately the coefficients of cos(4(—^)^/^/3), sin(4(—C)^/^/3) and the “non-oscillatory” terms. 
One finds that the expansion is of the form 


(5.55) 

m) = (-C)'/" + 


+(- 0 -''“s!'”sm(|(- 0 ’''‘)+ 4'’»=( 

where S^\ S'® 

and are asymptotic series valid for ( — oo: 

(5.56) 

So°'’ = ^(^1 +wl + wl+ wi), 

(5.57) 

5'^ = ^{wl-wl-wl + wl) , 

(5.58) 

^( 0 ) _ (_C)- {w^W/g - W 1 W 2 ) ■ 


Now using the definitions ( 5.53 )-( 5.54 ) of the asymptotic series (with ( replaced by (—we see 

that is an asymptotic expansion in powers of (—C)~^j starting with the constant term. Similarly, 
is an asymptotic expansion in powers of (—C)~^ starting with (—and is an asymptotic expansion 
in powers o f th e form (—p = 0,1,.... This establishes the form of the expansion specified in Step 
1, formula (5.S) 

Proof of Steps 2 and 4 : No w we must establish the existence of and Gt, ^ > 1, satisfying (^.9|)-( ^T^ ) 
(for £=1,2, 3), and (|^ ) -(|^ ) (for £ > 4). 

We observe that if {Fi}£>i satisfies (^.15)-(^T7) and (5.32), then clearly Fg is the £-th antiderivative of the 


function Fq. Now Fo(C) = (^^^(0) ~ C (^*(0) > and it straightforward to prove that the £-th antiderivative 
of Fq is given by 

(5.59) FgiO = {At\0f+q[^\C) {M<:)f + qi^\OAimi'iO 
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where are polynomials in C. Now ( 5.59 ), togeth er wi th the asymptotics for ^ > oo of the Air y fun ction 

Ai{() imply that for each £, Ft satisfies ( 5.18| ) and ( ^.36| ). Furthermore, using the asymptotics ( ^.50| ) and 
( ^.51 ) (again replacing z with —C, and C with (2/3)(—and letting C ^ —oo, we find that for each £, 
Ft possesses an expansion of the following form 

(5.60) FtiO = (-C)'/2 + 


M-CF‘ d'> sin + sf SOS , 


where So, Si, and S 2 are asymptotic series valid for C ^ — 00 . We will refine our knowledge of these series 
in what follows, but for now we observe the following important basic property of these asymptotic series: 
Property 5.61: The asymptotic series Sq^\ s[^'^ and S^^'^ are all in integer powers of (—C)- 


In order to obtain more detailed information about the asymptotic series 5”^, S'C and S '2 j we will make 
use of the following Lemma. 

Lemma 5.1. If a function G(f) is C°° on (— 00 , a) for some a > 0, and possesses the following asymptotic 
expansion (which may be differentiated repeatedly to obtain asymptotic expansions for derivatives), 


(5.62) 


G(C) = (a(-C)^ + K-C)^ + ... ) tr^g (^(-C)"/") , 


in which the asymptotic expansion in parentheses descends in powers of (—j ^ tfl-gi') denotes 

either sinf) or cos (•), then the function G(()d( possesses the following asymptotic expansion for s ^ 
— 00 : 

GiOdC = + bi-s)^ + . ..)trtg (^(-5)'^^) 

where the expression in parentheses descends in powers of (—s)^/^. 

Proof. 

-(j + 3) 


(a(-C) 

(a(-C) 


-C “ + 


5(_0^ + ...)tr*5(^(-C)^/^) dC 


-O^+bi-0' 


-0 + 3) 


-1 


/ / , ^ \ - (j + 3) 

- fa(-s) 2 + b{-s) 2 + 


2(-C)i/2 dC 
1 


trig ( g(-C)^^^ ) df 


I trig 




(a(-C)' 


(j+i) 


+ &(-C)' 


(a(-s) 


j + 1 / -(j + 3) 


2 (-s)i/2 

‘T+ + ...)'t+s(^i(-0+’^<lC 

5 ' 

-(j + 6) 


~ (j + 1 ) — ( j + 4 ) \ /^/ \‘^/9 

2 +b{-s) 2 _\trig[-{-sft‘^ 




-1 


X 


'2(-C)l/2 df 


trig dC„ 


and continuing in this way establishes the lemma. 
Let us write Fi in the following convoluted way: 

(5.63) 


□ 


m) = J_ {ms) - coi-sy/^) ds - ^coi-cf^. 


One sees that Fi so defined is an antiderivative of Fq immediately. To verify that Fi so defined is exponentially 
decaying for f + 00 , one uses the asymptotic expansion of Fp obtained in Step 1 together with Lemma 5.1 


to compute the asymptotic expansion of the integral in (5.63), and compares to the asymptotic expansion 
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( ^.60D. S ince there is no constant term in (^.60 ), neither in the asymptotic expansion of ( [S.63| ), Fi defined 
by ( ^.63 ) must be expressible in the form (5.5£), and hence it is exponentially decaying for ^ > +oo. 

So we have learned that (1) Fi exists satisfying F[ = Fq, (2) Fi satisfies (5.18) for ^ ^ +oo, and (3) 
Fi possesses an asymptotic expansion of the form ( ^.60 ) for (( —oo. Furthermore, we learn that the 
asymptotic series (—appearing in ( 5.6(| ) with ^ = 1 is obtained by computing the term-by-term 
anti-derivative of the asymptotic series (—C)^^^'5'o appearing in ( [3.55 ). This provides us with more detailed 
information about the asymptotic series S'q^^: 


(5.64) 




If, rather than using Lemma 5.1, one is more careful integrating by parts repeatedly, it is straightforward to 
prove the following two properties: 


Property 5.65: The expansion 5'^^^ starts at (—()) ^ and descends in powers of (—C)^- 

Property 5.66: The expansion S' 2 ^^ starts at (—C)~^ and descends in powers of (—C)^- 

Setting Gi = Fi -I- (2co/3)(—we have proven the existence of Gi satisfying (5.£) and ( 5.12| ), with 

= — 2ci/3 (the constant Ci appears in ([3. 8 ])). We may now define 


(5.67) 


Go = 


Gi{s)ds 


The function G 2 clearly satisfies (5. H), and using the asymptotic expansion ( |5.12D satisfied by Gi together 
with careful integration by parts shows that G 2 satisfies the asymptotic expansion ( [3.13| ) for ( —> — 00 . We 
next define F 2 via 


(5.68) 


F2 = G2 + 


which, it turns out, is equivalent to the definition F 2 = Fi(s)ds. The function F 2 then clearly satisfies 
(|t|) and ( ^ 11 ). 

The functions G 3 and F 3 are defined in similar fashion. First we define 


(5.69) 


G3= f (^G2 {s) - c^o\-s) 

J —00 ^ ' 


Second, using the asymptotic expansion ( 5.13| ) satisfied by G 2 and careful integration by parts in the manner 
used to prove Lemma we deduce that G 3 possesses the asymptotic expansion (5.14). Third, we define 
F 3 in the obvious way: 


(5.70) 




Uniqueness of asymptotic expansions then implies that F 3 = F 2 {s)ds, and hence F 3 satisfies ([3.17|) and 

(Pi. 


The inductive argument goes as follows: suppose we have the existence of and Ge for i = 1,..., 3k, 
satisfying (5.9)-(5.18), and possibly (5.31) - (5.37) with j = k — 1 (if A: > 2). We define G 3 k+i and G 3 k +2 as 
follows: 


(5.71) 

(5.72) 


G3fe+i(C) = (Gskis) - ds 

J — 00 ^ ' 


G-ij+i - 


G3fc+i(s)ds. 


The asymptotic expansion satisfied by G 3 h (either ( 5.14 ) if /c = 1, or ( |5.33| ) with j = fc — 1 if A: > 2), together 
with a careful integration by parts argument as outlined in the proof of Lemma |5.l| shows that Gsfc-i-i satisfies 
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the following asymptotic expansion: 

(5.73) G3fe+i = (-C)'/" (-C)-' + • • •) 

+Gf+') sin cos (-C)'/"), 


where G^ and are asymptotic expansions: 


(5.74) 


(5.75) 


For 3fc + 1 even: 


G^/c+l) ^ 


N-(3fe+6)/2 


^(3fc+l) _ ^—(3/c+3)/2 ^ j{Zk+l) y.(3fc+l) ^+ * • • ^ 

For 3k 1 odd: 

G^/c+l) ^ ^_^^-(3fe+3)/2 ^^(3fc+l) ^pfc+l)(_^)-3 _^ . . .^ 

^(3fe+l) _ ^—(3/c+6)/2 ^j(3fc+l) y(3fe+l)^_C)~^ ' ' ') 

Similarly, G 3 fc +2 satisfies the following asymptotic description for C ^ — oo: 


(5.76) 


G3fe+2 = (-C)'/" (c(""+")(-C)-' +cf+") (-C)"" + • • •) 

+Gr«> Bi„ (i + Cg‘«> cos (I (-0='"). 


where Gg and are asymptotic expansions: 


(5.77) 


(5.78) 


For 3k + 2 even: 


G0/C+2) ^ ^_^^-(3/c+ 7)/2 j'^(3fc+2) ^(3k+2) ^ ^ ^ ^ 

^(3fc+2) _ ^(3/c+4)/2 ^^(3/c+4) ^ y.(3fc+4) ^+ * * * ^ 

For 3k 2 odd: 

G^3fe+2) ^ ^_^^-(3fe+4)/2 ^^(3fc+2) ^(3/c+2) ^ ^ ^ ^ 

^(3A:+2) _ ^—(3/c+7)/2 ^y.(3fc+2) ^ ^(3A:+2)^_C)~^ ‘ ‘ 


We may now define ^ 3 ^+ 3 : 

(5.79) G3fe+3(C)=r (G3fc+2(s)-c(""+")(-s)-G2)ds_2c('"+"^(-C) 

J — 00 ^ ' 


1/2 


It is by now a straightforward exercise to establish that G 3 k +3 so defined satisfies (5.31) as well as the 
asymptotic description (5.33) for C, —oo. 

Turning now to the existence of F^k+i, 7 ^ 3 /c +2 and F^k+s, we define them through 


(5.80) 


7 ^ 3 fc+l = / F 3 kis)ds, F3k+2 = / F 3 fc+i(s)ds, ^ 3^+3 = / F 3 k+ 2 is)ds. 

J 00 j OC j OO 


It is immediately clear from these definitions that {^ 3 ^+^}^^ satisfy (5.18) for ( —> 00 . 
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Because of the recursive definitions of Fi and G(, it follows that 


(5.81) 


Fsk+i — Gsk+i — 


+ - 


2 c 


{3k) 


3/2 


3 

(3fc-3) 


945 


(-0 

(- 0 ='= + 


+ 


(_1^3fc+l23fc+l^^ ^_^^(6fe+3)/2 


n-=l'(2j + i) 


(5.82) 


^3fe+2 — G3k+2 + 
2 ^ 


4c 


(3fe) 


15 

j3fc-3) 
'0 


(-0 


5/2 


(5.83) 


10395 

d3fe) 


(-0 


11/2 


(_ l ) 3 fc + 223 fc + 2 ^^ ^ 06 fc + 5)/2 


n-=f(2j + i) 


(- 0 ^ 


F. 


3fe+3 


= G; 


8 c 


3A;+3 


0 


+ - 


2 ®, 


105 

,(3fc-3) 


(-0 


7/2 


135135 


(-0 


13/2 


(_ l ) 3 fc + 323 fc + 3 ^^ (_^)( 6 fc + 7)/2 


n?=f( 2 j + i) 


we observe that there are no constant terms in the relations ( ^.81 )-( ^^ ). Indeed, this is so because (1) for 
each £, Fg has the representation ( 5.59| ) for some polynomials and consequently the asymptotic 

no constant term and ( 2 ) Gi as defined 
-*■— 00 . Therefore the arbitrary constants 


expansion for F^ valid for C —> — oo is of the form (5.60) with 
possesses an asymptotic expansion with n o con stant term for - 
of integrati on w hi ch m ay appear in ( ^.81 )- (|5.83| ) are all 0. _ 

Finally, (^.8l|) - (^.83| ) imply that {T 3 fc+M}^=i and {G 2 k+ti}%^i satisfy (F^ ) with j = k. 

This completes the proof of Step 2 an d Ste p 4, and we have established the first claim of Theorem 
that asymptotic expansions of the form ( |l.35| ) hold true. The explanation for why the coefficients in such 
an asymptotic expansion are analytic is as follows. 

We return to the integral (5.4), and show that each term in its asymptotic expansion depends analytically 
on t. Since the same is true of the other integral appearing in (5.3), this will complete the proof of Theorem 


1.3. We begin by observing that each term in the asymptotic expansion of (5.4) is of the form 

r0 


xfiX)Q{X,t) (-0/3)'/'dA, 


where 


• Xp^X) represents some finite number of derivatives of X/ 3 (A); 

• Q{X, t) represents terms obtained from the integration by parts procedure but do not depend on the 
partition of unity X/S (these terms depend analytically on A and on t); 

• may be taken to be analytic with branch cut emanating from X = (3, down the real 
axis, and passing in particular through z* + e. For A away from the cut, this function also depends 
analytically on t. 

We now decompose the integral into a sum: 

r0 


1/2 


dX. 


(5.84) r (A)g(A,t)(-0/3)'/'dA= / x?(A)Q(A,t)(-0/3)'/'dA+ / Q(A,t)(-0/3) 

The first integral on the right hand side of ( ^84 ) is clearly analytic in t, because the integrand is analytic in 
t. For the second integral, note that X /3 does not appear. This is because the partition of unity is such that 
X/3 = 1 on {z* + e,/3). Now this second integral is seen to be analytic by expressing it as 1/2 the value of 
a contour integral encircling the interval {z* + e,/3), passing through z* + e. This can be done because the 


integrand is analytic, with a square-root branch point at A = /3. This completes the proof of Theorem 1.3 











































42 


N. M. ERCOLANI AND K. D. T-R MCLAUGHLIN 


6. Appendix 


Here we present a set of explicit formulae for S'i(z), the solution of Riemann-Hilbert problem 3.58. 
1 . ioT z€C\{BfUB^) 

(6.1) S,iz) = 


144/i(a) 


7 


48(a — (3)h{a) 


3(2-a) " + - 
{z-a)-^ 


1 _ {a) 


+ 


144h(/3) 

7 


3(2-/?)-" + - 


a — f3 h{a) 
— 1 i 
i 1 

1 


{z-a) 


-1 


-3 


h'ip) 


48(/3 — a)h{(3) 


(z-p) 


-1 


5 iP — a h{p) 
— 1 i 

i 1 




— 1 —i 
—i 1 


-1 i 
i 1 


2. For 2 G , 


( 6 . 2 ) 


-S'l( 2 ) = 

5 

144h(a) 


7 


48(a — P)h(a) 


3(2 - a) ^ + 
( 2 - a)- 


1 _ 3 ^'(a) 


72 i2h(/3) 




a — P h{a) 
— 1 z 
i 1 

3 


( 2 -a)- 


— 1 —i 
—i 1 


10h{P) 


1 3h'{P) 


P — a h{P) 
(2 — a)^/^ 


{z-py 


7 

+ 722 


72 2{P-a)h{P) 


(z-p)-^- 


(2 — /3)^/^ Jp R{s)h{s)ds 

{z-pyi^ 


— 1 i 
i 1 

f 1 i 


{z — ay/‘^J^R{s)h{s)dsj V * 


3. For 2 G 

(6.3) ^ 1 ( 2 ) = 

5 


144/i(/3) 


3(2-/3)-2 + ^ 


1 


7 


48(/3 — a)h{P) 


{z-py 


72 {2h{a) 


{z — a) ^ + 


P — a h{(3) 
1 z 
z —1 

3 


{z-py 


— 1 i 
i 1 


1 3h'{a) 


lOh(a) [a — P h{a) 

{z-py/^ 


{z-a) 


-1 


-1 -i 


{z — R{s)h{s)ds \ V ^ 


^jy 


-{z-a) ^ - 


(2 — a)^/^ 


1 —z 


72 2(a -/3)/i(a) {z - P)^/^ R{s)h{s)ds j \ -i -1 

For comparative purposes, if tj = 0,j y 4, and if we set t 4 = t, then we have 
1. for2GC\(H^URf) 

(6.4) S'i( 2 ) = 

-7/3 /3(-562 - 136/3+ 15/3^2 + 25/3^) W 1 i 


(6.5) 

( 6 . 6 ) 


+ 


96(8-/32)(/3-2) 
7/3 

96(8 -/32)(z + /3) 


96(8 -/32)2(z + /3)2 


z —1 


/3(-136/3 + 562 + 25/33 - 15/32z)W -1 / 


96(8 -/32)2(/3 -z)2 


z 1 
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(6.7) 

( 6 . 8 ) 

(6.9) 

( 6 . 10 ) 


-7(3 


\96{8 - P^){z - (3) 

/3(-56z - 136/3 + 15/32 z + 25(3^) 
96 ( 8 -/ 32)2 ( 2 + / 3)2 

7/3 


77 " 


72 R{s)h{s)ds ^ 


1 i 
i —1 


96(8 -/32)(z + /3) 

/3(-136/3 + 56/3z + 25/3^ - 15(3^z) 


3. For z G Bg we have 


( 6 . 11 ) 


S'i(2:) = 


— 1 / 

96(8 -/32)2(/3 -z)2 72-f^ R{s)h{s)ds J y i 1 

-/3(-136/3 + 56z + 25/33 - 15/3^2) 


9Q(8- P‘^Y{I3- zf 


7_ 


47 


-2 


-1 


_ I3{z + I3) 

3 R{s)h{s)ds (8 — /32) 


— 1 / 
/ 1 


(3 (25/33 - 136/3+15/3^2-56 z) 
96 (/32 - 8)^ (z + /3)^ 


72 J^gR{s)h{s)ds 


7_l 

96 


8-/32 


-1 -/ 

—/ 1 


7. Conclusions 


This paper provides the first mathematical derivation of the form of the asymptotic expansion of the 
Hermitian random matrix partition function, conjectured in |^. This uses the connection to non-classical 
orthogonal polynomials, together with the Deift-Zhou steepest descent and stationary phase method for the 
asymptotic analysis of Riemann-Hilbert problems. In particular, the asymptotic analysis of 0 presented 
here demonstrates that 

• log ( Zn ) has an asymptotic expansion of the form of 


(^); 


• For arbitrary integers g > 0, eg(ti,... ,t,^) is an analytic function of the (complex) vector t := (3i,... ,t^), 
in a neighborhood of ( 0 ,..., 0 ); 

• For times, t, within the domains of holomorphy for the eg, the limiting mean density of eigenvalues for 

the Hermitian random matrix ensemble is supported on a single interval {a, (3). The computed asymptotic 
expansion of valid for all A G (a, (3) and in fact depends only on the equilibrium measure; 

• the coefficients of the Cg are generating functions for graphical enumeration in the sense that 


d 

= #{k-valent, n-vertex, g-maps}. 


As a consequence of these results one is now situated to be able to combine asymptotic analytical tech¬ 
niques with the results of many ingenious ideas, contributions, and calculations which have arisen over the 
past 20 years from the theory of 2D quantum gravity 113 , to investigate the connection between the asymp¬ 
totic expansion of the partition function and solutions of integrable systems such as the KdV hierarchy, the 
Toda lattice, and singular limits of these equations described by modulation equations 1^, 1^ . 
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